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Chapter 1 

SOME GENERAL FACTS 



After a short introduction, the most important known facts about sine-Gordon and massive 
Thirring models are exposed. It is also explained their connection with the c = 1 free boson. 

1.1 Introduction 

As a very large number of papers in literature, this thesis principally deals with the sine-Gordon 
model, which has been well known at the classical level for the late fifty years and plays also an 
important role in quantum theory, thanks to its particular properties of non-linearity and inte- 
grability. It has been successfully applied in very different sectors of Mathematics and Physics, 
from partial differential equation theory to particle physics or solid state physics. Recent applica- 
tions of the classical model are related to nonlinear optics (resonant dielectric media) and optical 
fibers, magnetic properties of polymers, propagation of waves in crystals and so on. Interesting 
applications of the quantum model are related to Kondo effect and to the thermodynamics of 
some chemical compound, as can be found in |^ (see also section |4.10 ). At the same time, the 



quantum theory shows a phenomenology that is similar to the Skyrme model used before QCD 
era to describe barions and strong interactions. 
The most relevant properties of the model are 

• at a classical level, all the solutions of the equations of motion are known (exact integrability 
via inverse scattering method) 

• the classical solutions describe solitons, antisolitons and bound states (breathers)Q ; in a 
scattering process this solutions are transparent (it is the mathematical meaning of "soliton") 

• it admits, both at a classical and at the quantum level, a countable infinite set of conserved 
charges 

• the quantization of the theory describes an interacting particle with its antiparticle and, in 
a certain (attractive) regime, bound states 

• the S matrix has been exactly determined; only elastic scattering processes can take place 
(i.e. no particle production), that is the quantum analog of the classical transparency of 
solitons 



^ The so called mesonic solutions are excluded in this analysis. 
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Finite size effects 



Finite size effects are widely recognized to play an important role in modern statistical mechanics 
and quantum field theory. From a statistical point of view, it is known that no phase transitions 
take place in a finite volume system. For example, specific heat c(r), that is divergent at the 
critical point, if the system has finite size looses it divergence; one observes only a rounded peak, 
in the plot c(T) versus T. Moreover, there is only an interval around critical temperature Tc where 
the finite size effects are relevant. Out of this interval, they are negligible (because only near Tc 
the correlation length can be comparable with the size of the system). The interesting fact is that 
specific heat (and other critical quantities) have a scaling behaviour (i.e. varying the size L) that 
is fixed by the (infinite size) critical exponents (see |jl2|). This is a general fact: as argued in p^ , 
the UV behaviour of the scaling functions (see later) is fixed by the conformal dimensions of the 
operators that belong to the universality class of the critical point (i.e. the CFT describing the 
critical point of the statistical system). 

Also in quantum field theory interesting phenomena appear. If the space-time geometry is a 
cylinder of circumference L, Casimir effects change the energy of a two body interaction, because 



particles interact in the two possible directions, as shown in figure |LlJ (looking forward one can 
see his own back). Also new radiative corrections to a propagating particle may appear because 
of the closed geometry.A system on a finite volume has discrete energy and momentum spectra. 





Figure 1.1: Finite size effects 



Liischer ||13|, |lj| has show that the corrections to the free system eigenvalues of E, P depends 



on the scattering amplitudes (i.e. on the S matrix). This general fact can be used to extract 
informations on an infinite space QFT from a numerical simulation, that obviously is affected by 
finite size effects. 

Notice that a cylinder geometry can also be interpreted as a finite temperature field theory 
(temperature is L^^). Such a theory is related to the phase transition between confined and 
deconfined phase in QCD. 



1.1. INTRODUCTION 
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In general, a physical quantity (for example the free energy) that depends on the renormal- 
ization scale L (it parameterizes the renormalization flow) is a scaling function. The energy 
and momentum computed in ( p.58| , 3 .591 ) are scaling functions, because in the NLIE the scale L 
appears. 

There are many methods of investigation of the finite size effects. One of them, that has been 
very fruitful even for theories which are not integrable, is the truncated conformal space approach 
(see 1 15 1 and section |l]J), which is an intrinsically non-perturbative approximation method. It has 
problems of principal nature, coming from the fact that one does not have an analytic control of 
the spectrum, and of practical nature, because to reach a certain precision in the resulting energy 
levels one has sometimes to resort to very high truncation levels and introduce enormous matrices 
to diagonalize. 

For integrable QFT, there also exist exact analytic methods to compute the finite size effects, 
like e.g. the Thermodynamic Bethe Ansatz (TBA), which was used to calculate the vacuum 
(Casimir) energy The method was later extended to include ground states of charged sectors 
[p^ ]. More recently, using analytic properties of the TBA equations extended for complex values 
of the volume parameter, an approach to get excited states was proposed in [^. Their method 
to get excited states sheds light on the analytic structure of the dependence of scaling functions 
on the spatial volume and up to now was the only method developed to deal with excited states 
in perturbations of minimal models. Its main drawback is that to obtain the equation for a given 
excited state one has to do analytic continuation for each case separately, and a major part of this 
continuation can only be carried out numerically. Because of the complications of the analytic 
continuation, this method is limited at present to simple cases of integrable perturbations of 
Virasoro minimal models and some other perturbed conformal field theories. Similar results were 
obtained in |jl9|, |4ll] . 

This thesis reports on a novel approach to the excited states of IQFTs in finite volume, based 
on the nonlinear integral equation (NLIE) method, which has its origin in the so-called light-cone 
lattice Bethe Ansatz approach to regularize integrable QFTs. It was argued in that sine- 
Gordon theory can be regularized using an inhomogeneous 6-vertex model (or equivalently, an 
inhomogeneous XX Z chain). The NLIE was originally developed in this framework to describe 
the ground state scaling function (Casimir energy) in sine-Gordon theory in and it was shown 
that in the ultraviolet limit it reproduces the correct value of the central charge c = 1. Similar 
methods were independently introduced in Condensed Matter Physics by other authors |17]. 

The NLIE was first extended to excited states in |32| where the spectrum of states containing 
only solitons (and no antisolitons/breathers) has been described. Using an idea by Zamolodchikov 
[^] they also showed that a twisted version of the equation was able to describe ground states 
of unitary Virasoro minimal models perturbed by the operator *l'(i,3)- A framework for generic 
excited states of even topological charge in sine-Gordon theory was outlined by Destri and de 
Vega in However, there has been a contradiction between the results of the two papers, which 
was resolved in ||2^, |9|] where it was shown that it was related to the locality and the operator 
content of limiting ultraviolet conformal field theory (GET). Besides that, strong evidence for the 
correctness of the predicted spectrum was given by comparing it to predictions coming from the 
truncated conformal space (TGS) method, pioneered by Yurov and Zamolodchikov in |15| and 
extended to c = 1 theories in |27, ^. Later a modification of the NLIE to describe the states of 
sine- Gordon/massive Thirring theory with odd topological charge was conjectured |28|. 



The NLIE for sine-Gordon theory was generalized to models built on general simply-laced 
algebras of ADE type in |42| for the case of the vacuum. More recently, in |^] P. Zinn- Justin 
extended the method to the spectrum of excited states for these models and he also made a first 
attempt to describe perturbations of minimal models of GET. A general framework for describing 
general excited states of minimal models perturbed by ^*(i,3) (only massive case) can be found 
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in ||3J], where is stated the correct form of NLIE equation to deal with this case, and also the 
simplest examples of the resulting excited states are checked (see chapter ^). 

In the following chapters, the general setup of NLIE will be presented, with the most relevant 
examples. 

Chapter 1 is devoted to summarize some well known facts that will be used in the following. 
In Chapter 2 the light-cone lattice is introduced and the Bethe equations for the 6 vertex 
model are written. 

Chapter 3 is devoted to obtain an integral equation equivalent to Bethe Ansatz, but that allows 
a continuum limit procedure. The analysis of the so obtained continuum theory is in Chapter 4. 

1.2 c=l CFT: free boson 

To fix some conventions and to define certain objects which are used later, a brief summary is 
given of the c = 1 free boson with a target space of a circle of radius R. The Lagrangian of this 
CFT is taken to be 



^ [ df,^ix, t)d''^ix, t)dx , X e [0, L] , (1.1) 



where L is the spatial volume (i.e. the theory is defined on a cylindrical spacetime with circumfer- 
ence L). In the sequel often the complex Euclidean coordinates will be used z = e27r(t-«x)/L^ ^ _ 
g27r(i-i-jx)/L rpj^g superselection sectors are classified by the C/(l)^ x ?7(1)^ Kac-Moody symmetry 
algebra, generated by the currents 



J{z) = idzip , J{z) = idzf ■ 
The left/right moving energy-momentum tensor is given by 



Tiz) = —dzifd^y, = V LkZ-''-^ , f{z) = —d,ipd,ip= V Lkz~ 

OVr ^ — ' OTT ^ — ' 

fe= — oo fe= — oo 



k-2 



The coefficients and of the Laurent expansion of these fields generate two mutually com- 
muting Virasoro algebras. If the (quasi)periodic boundary conditions are required 



(p{x + L, t) = ip{x, t) + 2TTmR , m G Z , 

then the sectors are labelled by a pair of numbers (n, m), where ^ (n is half integer because of 
the locality, see later) is the eigenvalue of the total field momentum ttq 

TTo = / vr(x, t)dx , tt{x, t) = —dtip{x, t) , 
Jo 47r 

and in is the winding number, i.e. the eigenvalue of the topological charge Q defined by 



2ttrI 



Q = t; — 75 / dx(p{x, t)dx . 



1.2. C=l CFT: FREE BOSON 
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In the sector with quantum numbers (n, m), the scalar field is expanded in modes as follows: 

(p{x,t)= (f){z) + ^{z), 

^z) = -ipo - ip- log z + i ^ k"'''^''' ' 

where the left and right moving field momenta p± (which are in fact the two U{1) Kac-Moody 
charges) are given by 



n 1 

P±=R^2'^^- (1-2) 



The Virasoro generators take the form 



~ 2 ' ■ ' ~ 2 ' ■ ' 

k=—oo k=—oo 

where the colons denote the usual normal ordering, according to which the oscillator with the 
larger index is put to the riglit. 

The ground states of the different sectors (n, m) are created from the vacuum by the (Kac- 
Moody) primary fields, which are vertex operators of the form 



V{n,m){z,z) =: e^pi{p+^{z) +p_^(2;)) : . (1.3) 

The left and right conformal weights of the field V(n,m) (i-e- the eigenvalues of Lq and Lq) are 
given by the formulae 

A± = f . (1.4) 

The Hilbert space of the theory is given by the direct sum of the Fock modules built over the 
states 



|n, m) = V^n,m)(fi,0) \vac) , (1.5) 
with the help of the creation operators a_fc , a_fc k > 0: 

n = ^ {a-fei ...a^kpa-h ■■■a-ijn, m), ki,... kp, h,... IqE Z+} 

(n, m) 

The boson Hamiltonian on the cylinder is expressed in terms of the Virasoro operators as 

HcFT = (^0 + ^0 - , (1-6) 
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where the central charge is c = 1. The generator of spatial translations is given by 



2tt 

P= — (Lo-Lo) . (1.7) 

The operator Lq — is the conformal spin which has eigenvalue nm on the primary field V(„^ • 
One can also introduce twisted sectors using the operator T that performs spatial translations 
by L: x ^ x + L. The primary fields V(^n,m) as defined above satisfy the periodicity condition 
TV(^n,m) = ^(n,m)- If the more general twisted boundary condition labelled by a real parameter v 
is required 

'TV{n,m) = exp {ivQ) V(„^^) , 

then it is possible to generate superselection sectors for which n E Z + 

It is important to stress that a particular c = 1 CFT is specified by giving the spectrum of the 
quantum numbers (n, m) (and the compactification radius R) such that the corresponding set of 
vertex operators (and their descendants) forms a closed and local operator algebra. The locality 
requirement is equivalent to the fact that the operator product expansions of any two such local 
operators is single valued in the complex plane of z. This condition, which is weaker than the 
modular invariance of the CFT, is the adequate one since the theory is considered on a space-time 
cylinder and do not wish to define it on higher genus surfaces. 



By this requirement of locality, it was proved in [20| that there are only two maximal local 



subalgebras of vertex operators: Ab generated by the vertex operators 

{^(n,m) : n, m G Z} , 

and Af generated by 

{y(n,m) ■ n e Z, m G 2Z orn G Z + ^, m G 2Z+ 1}. 

Other sets of vertex operators can be built, but the product of two of them gives a nonlocal 
expression. 

1.3 Sine-Gordon/massive Thirring field theory 

The minkowskian lagrangian of sine-Gordon theory is given byQ 

CsG = I Qs.cDS'^cD + ^ : cos (/JcD) :^ dx , (1.8) 

where $ denotes a real scalar field, while that of the massive Thirring theory is of the following 
form: 

CmTh = j [^{ilvd'' + mo)^ - ^^-i'-^^-i,-^) dx , (1.9) 



Integration sign without any specification means integration on the whole real axis. 



1.3. SINE-GORDON/MASSIVE THIRRING FIELD THEORY 
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Figure 1.2: The family of vertex operators V(^n,m) with n £ Z/2 and m £ Z. Sector I is the common 
subspace. I and III are Af, that defines the UV behaviour of massive Thirring; I and II are Ab that is 
the UV of sine-Gordon, IV is a sector of non mutually local vertex operators. 



describing a current-current selfinteraction of a Dirac fermion It is known that the two theories 
are deeply related provided their coupling constants satisfy 

_ 1 

An 1 + g/n 

For comparison with the Destri-de Vega nonlinear integral equation, it is important to deal with 
cylindrical sine-Gordon and massive Thirring, i.e. the integrals in (L8, p!^ ) must be taken in the 
interval [0, L]. 

The cos term in (1.8) can be considered as a perturbation of the c = 1 free boson compactified 
on a cylinder, as described in section L2. Similarly the massless 
c 



(mo 



0) Thirring model is a 

1 conformal field theory and the mass term plays the role of a perturbation. From Coleman's 



paper it is known that correlation functions of the perturbing fields and : cos /?$ : are 
identical, then both models can be considered as the perturbations of a c = 1 compactified boson 
by a potential V : 



HsG/mTh = HcFT + ^ , V = X (V(i^o) (z, z) V(_i^o) z)) dx 

Jo 



(1.10) 



which is related to the bosonic lagrangian ( |L^ ) by the following redefinitions of the field and the 
parameters: 



/3 ' 2/32 



(/. = V47r$, R=^, A = ;^. (1.11) 



For later convenience, a new parameter p can be defined by 
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The point p = 1 (i.e. 5 = 0) is the free fermion point, corresponding to a massive Dirac free 
fermion. The particle spectrum of sG for p > 1 is composed by the sohton (s) and its antiparticle, 
the antisohton (s). It is known as repulsive regime because no bound states can take place, p < 1 
is the attractive regime, because s and s can form bound states that are known as breathers. The 

values p = — , k = 1,2, . . . are the thresholds where a new bound state appears. The potential 

k 

term becomes marginal when = Svr which corresponds to p = 00. The perturbation conserves 
the topological charge Q, which can be identified with the usual topological charge of the sG 
theory and with the fermion number of the mTh model. 

Mandelstam |^] showed that a fermion operator satisfying the massive Thirring equation of 
motion can be constructed as a nonlocal functional of a pseudoscalar field (boson) satisfying the 
sine-Gordon equation. But the fermion and the boson are not relatively local and then do not 
create the same particle (the two theories are not equivalent). 

The difference between them is that they correspond to the perturbation by the same operator 



of the two different local c=l CFTs Ab and Af as in[L^. The short distance behaviour of the sG 
theory is described by the local operator algebra At, while the primary fields of the UV limit of 
mTh theory are Af. 

Note that the two algebras share a common subspace with even values of the topological 
charge, generated by {V(^n,m) '■ n G 7,, m £ 2Z}, where the massive theories described by the 



lagrangians (L8) and ( |1.9| ) are identical. Exactly in this subspace holds the well known proof by 



Coleman about the equivalence of the two theories ||21||. The figure |L2| shows the four sectors 
where all the vertex operators live. 

1.4 Truncated Conformal space at c = 1 

In this section is given a brief description of Truncated Conformal Space (TCS) method for c = 1 
theories. It will be usefull, in chapter ^ to have a comparison of the data obtained from the 
non-linear integral equation that will be introduced. 

The TCS method was originally created to describe perturbations of Virasoro minimal models 



in finite spatial volume ||l^. Here is presented an extension of the method to study perturbations 
of a c = 1 compactified boson, more closely the perturbation corresponding to sine-Gordon theory. 

The Hilbert space of the c = 1 theory (on a cylinder) can be split into sectors labelled by 
the values of P and Q, which are quantised by integers. The numerical computations shall be 
restricted to the P = sector (it is not expected that any relevant new information would come 
from considering P / 0). The TCS method consists of retaining only those states in such a sector 
for which the eigenvalue of Hcft is less than a certain upper value Ecut, so the truncated space 
is defined as 



nTCs{s,m,E,^t) = {\^) : P\^) = s\^) ,Q\^) = m\^) ,HcFTm <Sc«t|^)} • (1.13) 

For a given value of s, m and Ecut this space is always finite dimensional. In this space, the 
Hamiltonian, represented on the basis of the eigenstates of the unperturbed c = 1 Hamiltonian 
(i.e. eigenstates of Lq + Lq)^, is a finite size matrix whose eigenvalues can be computed using a 
numerical diagonalization method. The explicit form of this matrix is the following: 
^ In this section, an hat on a letter means matrix 



1.4. TRUNCATED CONFORMAL SPACE AT C = 1 
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where Lq and Lq are diagonal matrices with their diagonal elements being the left and right 
conformal weights, / is the identity matrix, 

, = AJ + Ap = £ = ^ (1.15) 

is the scaling dimension of the perturbing potential V and the matrix elements of B are 

B^,^ = \ (^1 V(i,o)(l, 1) + l^(o,i)(l, 1) 1^) • (1-16) 

The units are chosen in terms of the soliton mass M which is related to the coupling constant A 
by the mass gap formula^ obtained from TEA in [|25|| : 



\ = kM^~^, (1.17) 



where 



2r(/i/2) 

7rr(l - /i/2) 



2r 



fl.15 



In what follows the energy scale is normalized by taking M. = \\ the dimensionless volume M.L 
is denoted by I. For numerical computations, the dimensionless Hamiltonian 

will be used. The usefulness of the TCS method lies in the fact that it provides a nonperturbative 
method of numerically obtaining the spectrum (the mass gap, the mass ratios and the scattering 
amplitudes) of the theory. Therefore it can serve as a tool to check the exact results obtained 
for integrable field theories and get a picture of the physical behaviour even for the nonintegrable 
case. The systematic error introduced by the truncation procedure is called the truncation error; 
it increases with the volume L and can be made smaller by increasing the truncation level (at the 
price of increasing the size of the matrices, which is bound from above by machine memory and 
computation time). 

Let us make some general remarks on how the TCS method applies to c = 1 theories. First 
note that the Hilbert space (even after specifying the sector by the eigenvalues of P and Q) consists 
of infinitely many Verma modules labelled by the quantum number n. At any finite value of Ecut 
only finitely many of such Verma modules contribute, but their number increases with Ecut- As 
a result one has to deal with many more states than in the case of minimal models. The results 
of TCS are supposed to approach the exact results in the limit Ecut oo. The convergence 



Notice the analogy with both s-G/mTh that are obtained perturbing the free boson theory by a "mass term". 
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can be very slow, while the number of states rises faster than exponentially with the truncation 
level. The perturbing operator has scaling dimension which ranges between and 2, becoming 
more relevant in the attractive regime, while the number of states corresponding to a given value 
of Ecut becomes larger as moving towards p = 0, which affects the convergence just the other way 
around. 

Generally, the energy of any state goes with the volume L as 



M 



vr c 



12(A* + A^)) 



6^ 



fc=i 



(1.20) 



where Aip (Ai^) are the left (right) conformal dimensions of the state in the ultraviolet limit, B 
is the universal bulk energy constant (the vacuum energy density) and the infinite sum represents 
the perturbative contributions from the potential V. 



The bulk energy constant has also been predicted from TBA and reads |25| 



1 . /pvr\ 



^ = -4tan^£^j {1.21) 

(the same result was obtained from the NLIE approach in |^). This is a highly nonanalytic 
function of p and it becomes infinite at the points where p is an odd integer. In fact, at these 
points there is a value of k for which k{2 — h) = 1, and Ck (^') — > oo. The infinite parts of 
B and (^') exactly cancel, leaving a logarithmic (proportional to I logZ) and a finite linear 
contribution to the energy, by a sort of a resonance mechanism. All of these "logarithmic points" 
are in the repulsive regime. However, due to UV problems in the repulsive regime we are not able 
to check numerically the logarithmic corrections to the bulk energy. 

The origin of UV divergences can be understood from conformal perturbation theory (CPT). 
It is known that when the scaling dimension h of the perturbing potential exceeds 1, CPT suffers 
from ultraviolet divergences which should be removed by some renormalization procedure. The 
TCS method is something very similar to CPT: it operates in the basis of the UV wave functions 
as well, but computes the energy levels using the variational approach and therefore could be 
called conformal variation theory" (CVT). As a result, it is expected that there could be UV 
divergences for the range of couplings where h > 1 which is exactly the repulsive regime p > 1 
|26|. The numerical analysis has in fact shown that in the repulsive regime the TCS energy 
eigenvalues did not converge at all when increasing the truncation level. 

Fortunately, there exists a way out: since it is expected to find a sensible quantum field theory 
when the UV cutoff is removed, it should be the case that the relative energy levels £m{L) = 
E-qi{L) — Eyac{L) converge to some limit. This is exactly the behaviour that has been observed. 
Consequently, in the repulsive regime one can only trust the relative scaling functions produced by 
the TCS method, while in the attractive regime even the absolute energy values can be obtained 
(including the predicted bulk energy constant ( |1.21 ), which is completely analytic for p < 1 and 
thus logarithmic corrections are absent as well). 

Many numerical results show that the smaller the value of p is the faster the convergence 
is (with the understanding that in the repulsive regime by convergence of TCS we mean the 
convergence of the energies relative to the vacuum). On the other hand, even in the attractive 
regime the convergence is so slow that to get reliable results (which means errors of order 10~^ — 
10~^ for the volume / ranging from to somewhere between 5 and 10) requires to work with matrix 
of dimensions around 4000. This means that the TCS for c = 1 theories is far less convergent 
than the one for minimal models (in the original Lee- Yang example the authors of |15| took a 17 
dimensional Hilbert space (!) and arrived to very accurate results). 



Chapter 2 

LIGHT-CONE LATTICE QFT 



In this chapter, the most important tools to deal with Destri-de Vega equation are introduced, 
with reference to the original papers [§, |6|]. Particular attention will be taken to indicate the path 
that has been done from the light-cone lattice until the definition of a particular system (the 6 
vertex model, alias XXZ chain), whose Hilbert space is quite completely known. 

2.1 Kinematics on light-cone 

It is a usual way to regularize quantum field theories by defining them on a space-like "hamiltonian" 
lattice (where time is continuous and space discrete) or space and time-like "euclidean" lattice 
(when both space and time are discrete). In statistical mechanics this is not just a regularization 
method but can be a right microscopic way to describe physical systems. In two dimensions, the 
most known approach is to define a rectangular lattice with axis corresponding to space and time 
directions and associate to each site an interaction depending only on the nearest neighbouring 
sites. In this case the partition function can be expressed in terms of a transfer matrix. 

In what follows, a different approach is adopted: Minkowski and Euclidean space-time can, in 
fact, be discretized along light-cone directions. Light-cone coordinates are: 

x± = X ±t 

and the choice 

M = {x± = -^n±, n± € Z} 

defines a light-cone lattice of "events" as in figure ^Tl] (a) . They are spaced by a in the space and 
time directions and by a/\/2 in light-cone directions. At every event P & A4 there is associated 
a double light-cone (in the past and in the future) and only events within this light-cone can be 
causally connected (see fig. ^ (b)). Then, any rational and not greater than 1 speed is permitted 
for particles, in an infinite lattice. The shortest displacement of the particle (one lattice spacing) 
is realized at speed ±1 and corresponds, from the statistical point of view, to nearest neighbours 
interactions. Smaller speeds can be obtained with displacements longer than the fundamental 
plaquette, and correspond to high order neighbours interactions. In quantum field theory, these 
are nonlocal interactions. 

In the following, only the local case (nearest neighbours) is treated. The nearest neighbours of 
the event P are the four nearest points in the light-cone directions. This implies that particles can 
have only the speed of light ±1 and are massless. They are called right-movers (R) and left-movers 
(L). 
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(a) 



(b) 



Figure 2.1: (a) light-cone lattice; (b) double light-cone emanating from a point and high-order interactions 
(dashed) 

In this case, it is possible to introduce a useful language for connection with statistical me- 
chanics associating a particle to a link. Consider the two links in the future that come out from 
a event P. Particles R and L in P, by definition, are respectively associated to the right-oriented 
link and to the left-oriented link. In this way, the state of a link is defined to be the state of the 
point where it begins (also the opposite choice, of connecting a link with the site where it ends, 
can be done; it is simply a matter of convention). For example, if in a point O there is a particle 
R, one tells that the "right-oriented" link outcoming from it is occupied by R. This correspondence 
of points and links is possible because only local interactions are assumed, and it is useful because 
the counting of states is simpler. But the physically correct interpretation is that particles live on 
events, not on links. Links are the possible world lines for particles. 




Figure 2.2: Periodicity in space direction; in link language the lines of time must be intended as shifted 



up by a little bit, as in figure 2.3 



In the following, the lattice is assumed of a finite extent L = aN in space direction (A^ is the 
number of sites, counted in the space direction), with periodic boundary conditions, but infinite 
in time direction, as shown in figure In this way a cylinder topology is defined for space time. 
The Hilbert space of states in an event P is the tensor product 



of R and L space of states. The fact that particles can be classified in left and right does not 
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mean, in general, that the two dynamics are independent, as happens for example in conformal 
field theory. In what follows, exactly the interacting case will be treated. 

Call \aLi,aRi) the generic vector of a basis of TCi where i = 1,...,N labels the sites. The 
notation 



\a2i-i,a2i 



\oiLi, am 



is useful and not ambiguous (even number refers to right, odd number refers to left). The total 
Hilbert space is: 



TV 



W7V = (g)Wi 



4 = 1 

and a basic vector can be represented by 

|ai, "2) ® ••• ® \ct2N-l,a2N) = |ai, 02, a2Ar) S Wat. 

Note that in a sites lattice, due to light-cone, 2N labels are required. 

If at a given time t there is a line of sites, the particular characteristic of the light-cone is that 
at time t + a/2 there is another line of sites, but not equivalent to the previous one, because it 



is shifted. Only at t -|- a there is an equivalent line (see figure [2.21) . Then two different evolution 
operators can be defined, depending on the initial state: 



C/+ |ai, 02, a2N,t) = "2' ct2Af' * + "Z^) 
{/_ \a[, a'2, a'2j^,t + a/2) = \a'(, a2jy,t + a) 



(2.1) 



where the initial states are chosen as in figure Schrodinger form of equations of motion is 



U. 





s 






<«1 


'^2/^3 





t-i-a 

t+a/2 

t 



Figure 2.3: Partial evolution operators. The convention used here for the labeling of states will be used 
also in the following. 

used, for a state in Hilbert space.Q The global time operator can be chosen as 

U = U+U- or U' = U-U+ 

depending on the initial state. For a consistent quantum theory, both this operators must be 
unitary. This can be guaranteed if the assumption C/+f/+ = u'LU- = 1 is made, that is the 
elementar operators themselves must be unitary. 



^ In Schrodinger form, if |q) is a state vector, its time evolution is given by \ci,t) — f/ja,0) where U satisfies 



-i / dtH 

motion's equations: U — Te J (Dyson's series) 
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Another operator plays an important role and is defined as follows (the states are at a certain 
fixed time): 



(2.2) 



it corresponds to an half-space shift in the space direction, with exchange of right and left states(see 
the figure 0). Two applications of V give a shift by an entire lattice spacing, then V"^ is the lattice 




.2r4-i 



Figure 2.4: Half-shift operator 

space evolution operator. It is possible to write an expression for V in terms of the permutation 
operator: 

Pnm {\an) ® bm)) = \am) O jOn) , l^n) O \am) & Un ® Ti-m- 

that looks like: 

V = Pl2P23-P2N-l,2N (2.3) 

The verification is direct: 

Pl2P23---P2N-l,2N \oil,a2, ...,a2N) = ^2-^23 • • -P2Af-2, 2Af-l ■■■ , Ol2N -2, a2N , Ol2N -l) = 



Pl2P23---P2N-3,2N-2 Oi2N , a2N-2, Ol2N-l) 



\a2N, Oil, a2Af-i) • 



Also y is a unitary operator. To show this, one can use the fact that the permutation is self- 
adjoint. A more complete list of known properties of P is: 



P = p- 



P\ 



P^ = 1 



(2.4) 



Then, the adjoint ofV is: 



= P2N-l,2N..-P23P^ 
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and VV^ = 1 simply follows. The operators defined up to now have the following commutation 
rules: 



(2.5) 



In figure 2^, the first case is proved, by showing the equivalence of the two paths V'^U+ and 
U+V'^. "Mutatis mutandis", all the other cases can be simply obtained. Consequently, the two 
principal evolution operators, V'^ and [/, are commuting: 

[y2^ U] = [U, V^] = 0. 

As previously shown, they are also unitary. Thanks to all these properties, it is very natural to 
identify them as the exponential of the hamiltonian operator, and the exponential of the linear 
momentum: 



U = e 



-iaH 



y2 _ g-iaP 



(2.6) 
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U 



+ 



U+ t 



Figure 2.5: Commutation rule: ¥"^11+ — U+V"^ 



There are other two important operators, defined as: 



Ur = u+v 



(2.7) 



As clearly shown in figure 2.6 for one of them, they correspond to one step evolution in light-cone 
directions. They are commuting and give the expressions: 




V 



Ur=U+ V 



u. 



+ 



Figure 2.6: Right-oriented shift operator 



U = UrUl 
[Ur, Ul] = 



= UrUI 
UrUr = uIUl = 1 



then, using also {UM yields: 



UR = e 2. 



Ul = e 2, 



(2i 



2.2 Dynamics on light-cone 

A dynamics can be defined by giving all the amplitudes of the different processes that can take 
place. The fundamental assumption is that in every site a whole process can happen, in the sense 
that if \aLiOiR)^^ and \I3li I^r) out ^'^^ ^^e incoming and outgoing states in a certain site, they can 
be considered asymptotic states and the transition amplitude is an S-matrix element: 



out {Pr,(3l \aR,aL)in - S{aji,aL)„„{f3R,0L)out 



(2.9) 



This, in general, is an m ^ n scattering. The system is then defined via its microscopic ampli- 
tudes. 
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Prom the definition (p.ip of evolution operators, the following expression holds (the "in" state 



' out 



is at time t, the "out" state at t + a/2, with reference at the figure 2^: 

{ai,a2, ...,a2N,'t\ U+ "2) «2Ar> =in ("1,02, •••,a2Af| "'l,"2)---'«2Af)c 

=in (ai,a2 |a'l,a2>o„i •••m (a2Ar-l,a2Ar 1"^-!' «2Af)o„t = (2.10) 

= TT 5^ 

ll«-l,...,Af (a2i-l,02i)in.{"2i-l'"2i)o"* 



The operator U- can be obtained from ( |2.5| ). In the last line, the product is on all the sites at a 
given time. 

It is interesting to observe that there is a powerful connection with statistical mechanics, at this 
point. Assume, for simplicity, that the lattice is euclidean and has a very large (thermodynamics) 
but finite extension in the "time" direction. In this way there is an initial time (with an initial 
state |a)j„) and a final time (with a final state |/3)o„t)- The total transition amplitude is given by 
the product on all sites of the amplitudes, and the sum on all the internal states (7', 7"): 



out VJ\ ^lin - 

int states site. 



The important fact is that if S^'^y/ is real and positive, it plays the role of a Boltzmann weight 
for a statistical system whose partition function is given by the previous expression, that was 
conveniently called Z. It depends only on the external states. The first appearance of this 
correspondence is in 

In both the statistical and the particle interpretation, the integrability (that will be always 
assumed in what follows) can be obtained by requiring that the scattering amplitude S^^y or the 
Boltzmann weight satisfy the Yang-Baxter equation (see ||2|, ^). In this case a general m —>■ n 
particles scattering can be factorised in the product of "elementary" 2^2 particles scattering. 
Then at every site there will be associated one of this "elementary" amplitudes: 

o«t(/3L,/3ij|a«,aL),„ = 5f«;^^ (2.11) 

and "in" and "out" states are two particle states (as in the previous equation). In statistical 
language, usually the Boltzmann weight is indicated with R and there is a little change of notation 
with respect to scattering case: 

t}I3r,Pl — oPr.Pl ("o 19^ 



(the lower indices are interchanged). A pictorial interpretation of that is in figure 2/7 



dc 




t * 



= s 



ab 



a b 

Figure 2.7: Pictorial interpretation of R and S matrix 

The Yang-Baxter equation takes the form: 

'-'12'-'13'-'23 — '-'23'-'13'-'12 



(2.13) 
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(with the indicated change it holds clearly also for R-matrix). 

Apparently, this "phenomenological" approach is quite unusual, because in traditional lattice 
quantum field theory at every site is associated one interaction potential (ex: 0^(i) is the potential 
on the site i), not a whole scattering process. This can appear as a sort of "macroscopic" approach, 
not based on fundamental interactions. But the properties of factorisable scattering must be 
taken into account. Factorization of generic amplitudes in 2 — > 2 particles amplitudes is a sort of 
quantum superposition principle and the remarkable fact is that between one scattering and the 



other, the particles are asymptotic ones, that means that they are free. For example in figure 



(1+1 QFT on the continuum) between the point 1 and the point 2 the motion is free. Every point 




Figure 2.8: Factorizable scattering. Between two interaction points the particles are free. 

contain all the interaction. This is what was assumed in the definition of the lattice. Then it is 
perfectly justified that every site is connected with a whole 2 — > 2 particles scattering process.^ 
The previous equation (p.lO|) can be written now in a more specific form: 



-'i'aia2 -"-030:4 ••• -"'02JV-102JV {■^■■'-^J 

The physical system dynamics is therefore defined by the assignment of an R matrix, and inte- 
grability is guaranteed. The consistency on quantum mechanical interpretation is obtained by 
requiring unitarity, hermitian analyticity and crossing symmetry for the S matrix obtained by 

It is important to obtain an explicit expression for Ul and Ur : 

(01,02, ■■■,a2N,t\ Ur |a'i,a2, ...,a'2N,t) = ^ {'^,t\ U+ |Y) (Yl ^ 

~f 

= T.^ZZ RllZ ■■■ (7i,-,72iV |a'2^,a;,...,4^_i,t) = (2.15) 

_ o"2]v"'l TjO'Wa p"2JV-2"2JV-l 

and, in a similar way, 

(ai,a2, 0!2N,t\ Ul |a'l,a2, Ci2Ni't) = Ralai Ratal ■■■ Ral%'^la2N (2.16) 



This expression is consistent with the figure p.6| . The important fact is that this operators can be 
expressed in terms of the transfer matrix of an inhomogeneous lattice model, that will be defined 
in the next section. 



^ In the next paragraphs, it will be explained that, for the particular case of the 6 vertex R matrix, a more 
traditional lattice QFT approach can be formulated, in terms of a fermionic field. 
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2.3 Euclidean transfer matrix 



Consider a two dimensional euclidean square lattice, with periodic boundary conditions, in both 
the directions. The links are the physical objects of the system. They can be in different states 
belonging to the vector spaces A and V. In principle, this vector states can be different, but in the 
following they will be identified. At every site can be associated a Boltzmann weight depending on 
the four links crossing at this site and on a spectral parameter A (see the figure ^l9| ). The simplest 



1 

M 



b =w(ablap, X) 



N 1 



(a) 



a 



(b) 



Figure 2.9: (a) periodic (toroidal) square lattice; (b) Boltzmann weight for a site; the labels are: 
a,b£A, a,/3£V 

case is to take the same Boltzmann weights in all the sites, but more general configurations are 
possible. In the following an inhomogeneity Aj will be assumed, where i is the column index (all 
the sites on a column have the same inhomogeneity), in the sense that the Boltzmann weight in 
the column i is taken to be: 



w{ab\a(3, A — Aj). 

Also for the boundary conditions it is possible to assume more general configurations than the 
simplest one (i.e. toroidal b.c). Assume that between the column N and the N+1 (that is 1) 
there is a nontrivial seam line, in such a way that the Boltzmann weights on the column N (with 
respect to the normal ones) are given by: 



w 



(^)(a6|a/3. A) = e*'^^u;(a5|a/3, A). 



(2.17) 



This choice is made because only the link b cross the seam line. This are called twisted boundary 
conditions (obviously the choice uj = reproduces the periodic case). 

The partition function for such a model is given by (see the figure p. 9]): 



Z= '[[w{ab\a(3, X- Xi)...w{ab\af3, X- XN-i)w^^\ab\aP, X- Xn) (2.18) 



link states rows 



(observe the twisted Boltzmann weight in the last column). It is simple to show that it can be 
expressed, thanks to periodicity, in terms of a row-to-row transfer matrix. To obtain that, it is 
convenient to define an operator t that acts on both the horizontal and the vertical spaces: 



t{X) : A ^ A 
U(A) = (a|t(A)|6) 



(2.19) 
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defines the horizontal action; tab depends on the horizontal links and acts on the vertical space as: 

(2.20) 



tabW : V ^ V 
(a| tabW 1/5) = [tabW]ai3 = w{ab\af3, A). 



The twist can be obtained with an appropriate operator that changes the horizontal states by a 
phase after the seam line: 



\aN+r, 



Au}A 



e Ai 



(2.21) 



where u; is a real parameter (twist), ^ is a self-adjoint operator and n is an horizontal index (in 
the vertical direction there is no twist). A must be chosen diagonal on the whole horizontal space 
A, therefore it acts as a phase. Observe that the lattice euclidean structure is a normal periodic 
one, because the twist acts only on the quantum mechanical structure. For a generic operator Bn 
acting on the site n the consistency with (^.21 ) requires that the twist acts as: 



Bn+u = e^'^B^e-^'^^. (2.22) 
This expression can be used for the determination of the tab operator acting on the last column: 



N {a\ tW(A) 16)^,,, =N {a\ tW(A)e-^ 



Aojb 



N{a\t(^)iX)\b), 



(2.23) 



where the "positions" of the states and operators are emphasized; in particular observe that only b 
crosses the seam line. Using the identification ( ^.201) the expression (|2.17 ) can be simply obtained. 
The row-to-row transfer matrix is an operator t^^\X,{\i},uj) that acts on the vertical space: 



N 



V 



i=l 



in the following way: 

t(^)(A,{Aj,u.) = j:rW 



ai 



ai,...ajv 



taiaii^ - Al) ta^asi^^ " A2) 



taNaii^ — AAr)e 



iujb 



(2.24) 



The Tab is called monodromy matrix. The expression for the partition function is: 

M 



Z = Tr 



y(JV) 



t(^)(A,{A,},^) 



Define now a matrix acting on AfS^ A given by: 



R^a'aW = [tab{X)U = wiab\aP, A). 



(2.25) 



(observe that in the operator in (|2.20| ) the Latin and Greek index belong to different spaces, while 
in R matrix they belong to the same space and this differentiation of notation is redundant). 
An obvious generalization is needed for the last column. This definition (|2.25|) shows that the R 
matrix is defined as a Boltzmann weight, as required in the section |2.2|. In the following will be 



explained the importance of (|2.25D for the integrability. 

It is possible to show, at this point, that the transfer matrix ( |2.24 ) can be used to express the 
operators defined in (|2.15 , 2.16D . Consider the following form for the inhomogeneity: 



(2.26) 
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where is a positive real number, and calculate a matrix element of the transfer matrix on vectors 
of y(2Af) (clearly, in this case, the periodicity requires an even number of horizontal sites): 

(ai,a2,---,a2Af| i('^)(e,{A.})| a[,a2, ...,a2Ar) = 

— ^aia2(0)aia'j ^a2a3(2©)Q,2a^ ••• ia2]vai(20)«2]va2JV ~ 



ai...a2N 



p"2°3 T?'*4°5 p"2]V"l 



"2]V"1 
'a2jv_ia2]V 



that is exactly the expression ( ^.16 ). In an operatorial form this means that: 

t(2^)(e,{Aa) = c/i. 



(2.27) 



A similar picture holds for the operator Ur, but to take into account this case an assumption on 
R is necessary. Assume that R matrix is hermitian analytic, that means: 



<(-A) 



that in S matrix language is S'^{s) = S{s*). This is a well known property of an S matrix and, 
as explained in the section the same requirement is necessary for R. Then, as for Ul, the 
following expression can be obtained: 



t(2A^)(-e,{Aa)t = [/«. 

Here the adjoint conjugation is only on the vertical space.^ It is given by: 



(2.28) 



U(A)1 



a/3 



'w{ab\(3a, A)*. 



Transfer matrix can express the evolution operators. 



As introduced in the section |2^, if the R matrix satisfies the Yang-Baxter equation (|2.13| ), 
then it is possible to show that the transfer matrix commute with itself at different values of 
the spectral parameter. In this case the system is integrable. The same conclusion holds also 
for the inhomogeneous and twisted models. Assume that the R matrix satisfies the Yang-Baxter 
equation, that in the explicit form is: 



bib2 I 



Using ( 2.25 ) it is possible to obtain another form, specific for the t operator: 



(2.29) 



(2.30) 



For consistency, the R matrix must be of the regular type, that means that at the origin it is the 
unit operator (or proportional to) on Afi^ A: 



This equation fixes a normalization for the R matrix. Then, to fit with the equations ( p.27 , 2.28) 
a normalization factor is required. It is fixed by the unitarity requirement. The ( 2.30| ) holds also 
if an identical shift is performed on both the spectral parameters: 



The definition of adjoint is tat \ct)* = {ce\ tl^t \P) 



2.4. 6 VERTEX MODEL: MAIN RESULTS 
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(the same for because only the difference appears in the R term. This shift is an inhomogeneity 
for the t operator. 

Consider now a group of matrices g £ G of dimension dim^ x dim^ and determinant one, 
such that they commute with R: 



[g(S)g,R]=0. 

Then define a transformed ("gauged" or twisted) t operator: 



(2.31) 



t 



{g)ab ~ gactcb- 



The twist introduced in (|2.23| ) is the special case of ^ = [/(l)*^™-^. It is simple to verify that also 
the gauged satisfies the same ( ^.301) . This is because of the commutation relation ( 2.31 ). The 
group ^ is a symmetry group for the Yang-Baxter equation and for the vertex model defined in 
this way. With some simple algebra, using many times eq. ( p.30 ), it is possible to obtain the form 
for the monodromy matrix, defined as in ( 2.2^ ) but with a possibly different gauge gi in every 
site of a row (all the column has the same gauge) and clearly an inhomogeneity. If the following 
synthetic notation is introduced, G = {gi, ...,g]\f), then the final result is: 

R{^ - \T^^\'d) Tk^\'&')\ = \t^^\^') rif ^(^9)1 -R(^9 - ^')- 



Taking the trace on the horizontal space A yields: 



(f)(^,{A.)),4^^«{A.)) 







(2.32) 



(observe that the gauge must be exactly the same). At the various values of t[f^(i9,{\)) 
describes an infinite family of conserved charges. Then the system is integrable. This holds in 
particular in the inhomogeneous and the twisted cases. 

The interesting fact, at this point, is that in some cases the transfer matrix can be exactly diag- 
onalized by Bethe Ansatz methods. This gives an exact expression for eigenstates and eigenvalues 
of the operator U. 



2.4 6 vertex model: main results 

The theory developed until now is general and not referred to a specific model. The simplest non 
trivial case to take into account in the previous framework is the 4 x 4 R matrix, corresponding 
to the choice A = V. As shown in ||3|, the most general solution is the so called 8 vertex model. 
This name means that only 8, between the 16 entries of the R matrix, are nonzero. A special case 
is the 6 vertex model, for which many results have been obtained in the light-cone description: 
this will be the principal object of this dissertation. The R matrix has the form (lower index are 
rows and upper index are columns) 



/a 



V 



a J 



(2.33) 



As Boltzmann weights, all this nonvanishing entries must be real and nonnegative. But Yang- 
Baxter equation is solved for all the values of this variables.^ As explained in ( ^.2[) , the scattering 
interpretation is possible only assuming the properties indicated in that section. 



Observe that the upper and lower entries are the same; this is the assumption of parity invariance. 
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There is a well known mapping between vertex models and spin chains (see |^), i.e. the 
transfer matrix is the exponential of the quantum hamiltonian of the chain: 



t 



(N) ^ -H 



In the case of 8 vertex model, the hamiltonian is the XYZ(l/2) chain, while in the special case 
of 6 vertex, is the XXZ(l/2) chain. In what follows, this identification can be useful to interpret 
some facts connected with Bethe Ansatz. The XXZ(l/2) chain hamiltonian is given by: 



N 



j+ii 



i=l 



and 7 is the anisotropy. The a are Pauli matrices. The total z-component of the spin will play 
an important role in Bethe Ansatz. 

This is a statistical approach, but it is possible to give a particle interpretation to the same R 



matrix ( |2.33| ) on the light-cone lattice. 

The simplest case of particles obeying Pauli exclusion principle (fermions) and without internal 
degrees of freedom (color number) is assumed. This means that in an event only one particle of 
type R and one L at most can take place. In other words, one link has two states: empty or 
occupied. At every point there are four links, that means 16 possible configurations associated to 
it. In terms of events, these are the possible configurations that connect a point with the nearest 
neighbours in the future. 

Assume now that only amplitudes that conserves the total number of particles (R+L) are 
nonvanishing. This reduces to 6 the permitted configurations, as it is shown in the figure ^.10 . 



This is simply the 6 vertex model whose R matrix is written in (2.33). The assumption of 









03 



CO, 



CO, 



CO. 



1 ^2 3 4 5 

Figure 2.10: 6 permitted amplitudes. Dots are particles. 



CO, 



integrability for this amplitudes gives the general six- vertex model. The requirement of symmetry 
under parity transformation implies that 103 = lo^ and = ujq. The convention adopted in the 
figure U and in ([2^ , |2^ , |2^ ) shows that 



UJ5 = iOQ. 



Now, this R matrix can be written in an operatorial form, by defining a lattice chiral fermion 
ipA,n, with A = R,L, and n labels the sites. The anticommutation rules are the canonical ones: 

{'4'A,n,tpB,m} = 0, ^tpA,n, '>pB,m] = ^AB^nm- (2.34) 

This fermion has some interesting properties, that are exposed in the paper and are sketched 
in the following list: 

1. the R matrix and all the other operators Uany can be written in an operatorial form in terms 
of the fermion; 
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2. the lattice hamiltonian, in the free case coi = 002 = b = 1, c = 0, can be exphcitly written; 
by this, the following dispersion relation can be obtained: E = ±k; this is the dispersion 
relation for a free massless particle; the unusual fact is that it is monotonous, then there is 
no doubling of fermions; this is a consequence of the nonlocality of the hamiltonian 

3. the lattice hamiltonian admits a continuum limit — > 00 and a — > but with L = Na 
fixed; the locality is recovered in this limit; the continuum equations of motion are those of 
the massive Thirring model. The space is compactified on a cylinder. 

In the continuum limit, the massive Thirring model emerges as the field theory characterising the 
scaling behavior of the dynamics on the lattice (remember that L is finite; the scaling behaviour 
is understood in terms of this L). 



2.5 6 vertex model: Bethe Ansatz 



The assumptions of unitarity and hermitian analyticity will be taken into account, for the R 
matrix, and this requires that the variables in ( p. 33 ) must have the specific form: 



a = a('!9, 7) = sinh(z7 — "&), b = 7) = sinh 1), 



c{'&,'^) = isin7 



(2.35) 



The transfer matrix defined by this R matrix ( |2.24 , 2.25| ) can be diagonalized with Bethe Ansatz 
method. In terms of the spin chain, this means that there are two operators, usually indicated by 
B['d) and C{^), whose expression is known, and there is a "reference state"|^ \ such that: 



C{'&)\n) = 



and 



B{{}{)...B{{}m)\^) 



(2.36) 



for appropriate values of i9j is an eigenstate of the transfer matrix. The "appropriate values" of -dj 
can be obtained as the solution of a set of M coupled nonlinear equations, that are called Bethe 
Ansatz equations. In general, because of ( p. 32 ), the transfer matrix contains all the conserved 
charges, in particular the hamiltonian. Then the Bethe Ansatz eigenstates are also eigenstates of 
the hamiltonian. The Hilbert space of the theory and the action of conserved charges on it are 
then perfectly know. 

All this computations for the 6 vertex model were obtained in js], |6|]; the final results are 
written here for the eigenvalues of the inhomogeneous and twisted transfer matrix: 



r(??, G, w) = e^"^ [a(7? - 9) a(i9 + 9)]^ JJ 



+e-- [5(^?_e)6(?? + 9)]^n*ii 
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^ The "reference state" at this point is only a mathematical object. Physically speaking, it corresponds to the 
ferromagnetic state with all the spins up. 
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and the values of '&j are defined by the set of coupled nonlinear equations called Bethe Ansatz 
equations: 
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+ ivr] 



(2.37) 



zvr 



where 2N is the length of the chain and the number of sites in a row of the light-cone lattice. 
The 'dj are called Bethe roots, and in principle, can take any complex value. But there is a 
periodicity in their imaginary part: 

^2 



7r . 

'&j + —I 

7 



then only a strip around the real axis must be taken into account for the Bethe roots: 

^2 „2-\ 



X I 



vr vr* 
'2^' 2^ 



(2.38) 



(2.39) 



Moreover, only the range 



< 7 < vr 



will be examined. 

Prom Bethe Ansatz it is known that the whole spectrum of the theory can be obtained using 
all the Bethe configurations having M < N and 'dj / i?^ for every j / k. In general for a state 
with M roots the third component of the spin of the chain is 



S = N-M, 



(2.40) 



because every operator B{'di) counts as —1 spin. 

This XXZ(l/2) chain has 2 states in every site, then 2^^ states. Then the energy spectrum is 
upper and lower bounded. Changing the sign of the hamiltonian gives another permitted physical 
system. This means that there are two possible vacua. The first one is the so called ferromagnetic 
ground state, corresponding to M = that is the reference state |r2) . It has spin S = N. The 
second one is the antiferromagnetic ground state, that can be obtained with M = N and all the 
roots -di real. It has spin 5 = 0. 

In what follows, only the antiferromagnetic ground state will be considered, because it has 
one important property: in the thermodynamic limit {N oo) it can be interpreted as a Dirac 
vacuum (a sea of particles created by B) and the excitations on this vacuum behave as particles. 

The energy E and momentum P of a state can be read out by the transfer matrix eigenvalues 
by use of the equations ( 2.27 , 2.2^ ). The final form is: 
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(2.41) 



Other integrals of motion can be obtained in a similar way by transfer matrix (see ( p.32|) ). Observe 
that the second term in the transfer matrix expression vanishes because b{0) = 0. 

For future analysis, it is more convenient to express the coupling constant 7 in terms of a 
different variable p: 



P 



vr 

7 



< p < 00. 



Then all the expressions will be in terms of p. 



Chapter 3 



A NONLINEAR EQUATION FOR 
BETHE ANSATZ 



In this chapter the fundamental nonlinear integral equation (on the lattice) for the Bethe Ansatz 
is derived. In literature it is known as Destri-de Vega equation. It was obtained first in ||^, ^ and, 
in the final form, in On this equation the continuum limit is performed. 



3.1 Counting function 

It is possible to write the Bethe equations (|2.37| ) in terms of a counting function, that will be 
called Zn{'&). To do that, one important preliminar definition is needed: 

sinh -^{iiTV + "d) 

(l^i'&.u) = i\og . ^' <P{-^,^) = 

smh ^(ivri/ - §) 

where the primary interest is in values u = 1/2, 1 and the oddness on the analyticity strip around 
the real axis defines a precise logarithmic branch: the so called fundamental branch. In the 
appendix ^ the complete structure of the cuts and all the relevant properties of this function are 
exposed. The counting function^ can be defined by 

M 

z^(i?) = Nicpi'd + e, -) + <p{^ - e, -)] -Y,^i^- ^k, 1) + 2u, (3.1) 

fe=i 

used to express the logarithm of the Bethe equations, one obtains simply: 

ZNi^j) = 2TTlj, GZ + ii^, <5 = (M)w2 = (iV-5)„od2G{0,l}. (3.2) 



The number Ij plays the role of a quantum number for the Bethe basic vectors ( 2.36 ) and it must 
be chosen depending on the value of 6. Notice that 6 and u; play a similar role, because both 
produce a shift in the quantum numbers Ij (if w is absorbed in the definition of Ij): in the first 
case the shift is exactly tt, in the second case it is a real (possibly complex) number. This means 
that the variable 6 can be absorbed in uj but the most convenient choice is to use them both. 
Observe that Bethe roots can be obtained as zeros of the equation: 

1 + (-l)^e^^^(''^) = (3.3) 



^ The name counting function will became clear in the next paragraphs, where the counting equation will be 
obtained. 
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3.2 Classification of Bethe roots 



From Bethe Ansatz it is known that a Bethe state ( 2.36| ) is uniquely characterized by the set 
of quantum numbers ^ , M < N that appear in (|3.2|) . Notice that M < N means 



5 > 0. The values of -dj to put in (|2^ ) can be obtained solving Bethe equations. It is also known 
that only states with 

are necessary to form a basis for the space of states. It is a sort of fermionic character for Bethe 
states [p^ . 

As usual in many cases, Bethe roots can be either real or in complex conjugate pairs (for large 
A^). In the specific case ( p. 37 ), there is another possibility, due to periodicit y (|2.3^ ): if a complex 



solution has imaginary part Qrwd = ^{p + I) it appears as a single (in ( ^.3?!) it produces the 
left hand side real, then its complex conjugate is not required). It is called self-conjugate root. 
Remember now that the maximal number of real roots (M = N) describes the antiferromagnetic 
ground state. From the point of view of the counting function, a more precise classification of 
roots is required: 



real roots; they are real solutions of ( |3.2D that appear in the vector ( 2.36| ); their number is 

Mr; 



holes; real solutions of ( ^ ) that do NOT appear in the vector ( ^.36 ); their number is Nh; 



• special roots/holes (special objects); they are real roots or holes where the derivative Z'j^{'dj) 
is negative^ ; their number is Ns; they must be counted both as "normal" and as "special" 
objects; 

• close pairs; complex conjugate solutions with imaginary part in the range < l^m-i?! < 
7rmin(l,p); their number is Mc; 

• wide roots in pairs; complex conjugate solutions with imaginary part in the range tt min(l, p) < 
|9m7?| < 7r£±i; 

• self-conjugate roots (wide roots appearing as single); Qrwd = vr^ii; their number is Msc- 

The total number of wide roots appearing in pairs or single is M^y. The following notation will 
be used (sometimes) for later convenience, to indicate the position of the solutions: hj for holes, 
yj for special objects, Cj for close roots, Wj for wide roots. 

Complex roots with imaginary part larger than the self-conjugates are not required because 
of the periodicity of Bethe equations. This classification is not at all academic; it will play an 
important role in the physical interpretation of the final equation. A graphical representation of 



the various types of solutions is given in figure 3.1. An important remark must be made: from the 



definition itself of Zn ( p.lj ) it is obvious that only for states without complex roots the fundamental 
strip for (p (see ( [B.3P ), that is the largest strip around the real axis without singularities, is the 
fundamental strip for Zj\[. In all the other cases the analyticity strip for Zn is narrower, and 
depends on the imaginary parts of the complex roots. 

An important property follows from this classification: the Zn function is real analytic if u! is 
a real number 

Zn (r) = {ZNm* (3.4) 



^ The characteristics of this type of solutions will be completely clarified in the next sections. 
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Figure 3.1: The different types of roots and holes and their position in the complex plane, /i denotes 
7rmin{l, p}. The upper line at ^ = + 1) is the self-conjugate one. 



3.3 Counting equation 

It is possible to obtain an equation relating the numbers of all the various types of solutions. 
The path is simple, and makes use of the asymptotic values obtained in ([B.5D to calculate the 
limits 1} — > ±00 in Zn (|3.1|) . Observe that the term — -dj, 1), in the case of the wide roots, 
takes contributions from the horizontal strips next to the fundamental, as in figure [B.l| . This 
contributions depend from M\yi (number of wide roots below the real axis) and M\yi- (number 
of wide roots over the real axis). The asymptotic limits are then: 



p — 1 



2uj 



Zn(+oo) =Ntt + tt- — -S + 27rsign(p - 1) Mw\ + 
p + l 

p — 1 

Zn{—oo) = —Ntt — TT — -j-j'S' — 27rsign(p — 1) Miy-^ + 2lj 



(3.5) 



Observe that the number of self-conjugate roots is: M^y-f — = Msc- On the lattice, all the 
values for N, S are permitted. But in view of the continuum limit, is expected to be larger 
than S: N ^ S (also the number of holes and complex roots is expected to be much smaller than 
N). Then 

Zn{+oo) > Zn{-oo). 

The function is globally increasing (but not necessarily monotonous: in some points it can be a 
decreasing function). Define now two variables C±) choosing k± in such a way that they take value 
in the indicated interval: 



C± = ±Zj\[{±00) + TtS — 2TTk±, —IT < C± ^ TT 

Using the previous expressions for Z^- the result is: 

C± = 27r 



S , u) 
± - + 



p + I TT 



1 S W 

2 ^ p+1 ^ 
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where the new symbol [xj is the integer part of x, that is the largest integer smaller or equal to 
X. In this way, the following expression holds: 



Zn{+00) = 27rlmax + TT + (+ 
Zn{-00) = 27rlmin - VT - C- 



(3.6) 



where Imax is the largest quantum number satisfying /, 

Z^i-oo) 



Zi\[{+oo) 

max < 7^ and Ijnin is the smallest 

zvr 



quantum number satisfying I„ 



> 



27r 



. This very precise definition is required to take into 



account for the special roots/holes that can appear in the tails (this can happen when < for 
very large values of 

The total number of real solutions is the sum of real roots and holes, and in terms of Imax, Imin 
it can be written as: 



Using the equations ( p. 40 



Mr + Nh = In 



+ 1 + 2Ns. 



p^ ) and the expression ( |3.5| ) the result is: 

Nh - 2Ns = 2S + Mc + 2e{p- 1) Mw+ 



1 S 
2^ p+l 



UJ 
TT 



1 s 

2^ p+1 



UJ 
TT 



(3.7) 



It is the so called lattice counting equation. Remember that S is a nonnegative integer. In the 
case of u) = 0, it turns out that Nh is even (remember that Mq is the number of close roots, and 
is even). 

The most important fact is that in this equation doesn't appear the number of real roots. This 
fact, together to what will be explained in the next paragraph, allows to consider the real roots as 
a sea of particles (Dirac vacuum) and all other types of solutions (holes, complex) as excitations 
on this sea. 



3.4 Non linear integral equation (I) 



In this section an equation generating Zn will be obtained. The counting function (|3.l|) can be 
written in the following way: 



0(^? + e,^) + 0(^?-e,i) 



Nh 
k=l 



Mc+Mw Mr+Nh 
k=l k=l 



(3.8) 



In this case collects close and wide roots, and Xk the real roots and holes. Now, it is convenient 
to deal with the derivative of this expression: 



i,'(^ + e,l) + 0'(^_e,i; 



k=l 



Mc+Mw Mr+Nh 
k=l k=l 



(3.9) 
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Both the previous equations hold for t9 G C. Let x be a real solution of the Bethe equation. 
Assume Z'^{x) ^ and define a complex neighbour of x by a small <C 1: fj, = x + i'. Clearly 
(— l)''e'^^^^^ = —1. Consider the expression: 

1 + (_i)'5e*^iv{i;+i') ~ 1 + (^_-^^-^s^iZt,{x)+iuZ'^ix) _ -i^z'^ix) (3.10) 
then the following identity holds: 



1 (-l)^e*^^(^'+^)iZ^(x + 



I' 



12 -X V 1 + (_i)<5e*2iv(£+'^) 



+ ... (3.11) 



(the dots are regular terms in ii — x). From the Cauchy theorem and from ( p. 11 ), given an analytic 
function /(x) on an appropriate strip containing the real axis, yields: 

' Jr.27rifi-x Jt^2-ki^^' i + (_i)5e*^iv(A«) ^ > 

where Tj. is a anti-clockwise curve encircling x and avoiding other singularities of the denominator, 
i.e. other Bethe solutions (real or complex). It is always possible to find such a closed curve, 
because Bethe solutions are finite in number. An equation like ( ^.12|) can be written for all the 
real roots of ( |3.2D , Xk, V/c. The derivative (/>'(??, 1) is an analytic function, if poles are avoided, 
and applying to that the expression ( 3.12| ), the last sum in ( |3.9| ) becomes: 



k=l k=l ^ ' 



(3.13) 



The sum on the contours was modified to a unique curve F encircling all the real solutions x^, 
and avoiding the complex Bethe solutions (this is possible because they are finite in number), as 
in the figure |3.2| . 



X 



X 



X 



T1+ ^ 

' ^< — X X e \ X e e >^ 

Tl- I ^ r 



X 



X 



Figure 3.2: Contours for the integration. The crosses are roots while the circles are holes. 
Clearly the F curve must be contained in the strip 

< r/+, r/_ < min{7r, vrp, \ '^m Ck\ V /c} 
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Without loss of generality, assume that r]^ = rj^ = r/, and deform T to the contour of the strip 
characterized by tj. The regions at ±00 do no contribute because of the vanishing of 0', then the 
integral can be performed only on the lines 12 = x ±irj, where x is rea]|^: 



dx 



27ri 



-x + iri, 1)- / ■ ^ — - + 



(3.14) 



The first term at the right can be written as: 



1 -|- (_l)5g*^jv(a;-i»?) 

-{-lYe-^^^'^'^-^'^'iiZ'j^ix - i7]) 

I -|- (^ — ly ^-iZn (x-iri) 



+ iZ'^{x - irf) 



(3.15) 



Putting this in ( 3.15| ) the integral contribution from Z'j^{x — irf) is independent from rj because 
no singularities are crossed if ry — > and at the infinity the integrand is zero. The following 
convenient redefinition can be introduced (see ( B.4D ): 



K{x) 



(/)'(X,1) 



sm ■ 



27r 
>"Tl 



2vr Ap + D,,,^:^ 
P + 1 



cos ■ 



27r 



In general, for every complex value of -d: 



(3.16) 



j dx [5{x)Z'j^{'d -x)+ K{'& - x)Z^(x)] = N 



0'(^ + e,^) + </)'(^?-e,i; 



/ 



Nh Mc+Mw 

+Y^4,'{{}-hk,i)- '^'(^-6,1)+ 

fe=l k=l 

dx . _ -(-l)'^e-*^'^(^-^'')iZ^(x - ir?) 



- / —ct)'{d-x + ir],l 



+ 



1 _|_ (^ — \YQ-iZi^{x-ivi) 

dx , . (-l)^e^^^(^'+^'?)iZ^(x + i7y) 

— [d-x- iri, 1) ^ ^ (_i)5e*z^(x+ir,) 



+ 



The second term on the left takes the form 

dxK{^-x,l)ZNix) = I dxK{x,l)Z'^{'d-x); 



+ 



(3.17) 



that can be obtained by shifting the integration line: x — > x + Qm^. But if poles are crossed 
(this can happen if the imaginary part is sufficiently large), their contribution must appear in the 
right hand side of the equation. 

^ f = In all the cases where the limits of integration are not explicitly specified, the integral is taken on 

the whole real axis. 
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Using Fourier transformations, the convolution in the left hand side of ( 3.17 ) can be simply 
handled. Observe first that the Fourier transform of [6{'& — x) + K{'& — x)] is simply 1 + K . In 
the appendix |B| the form of K is given (as cj)' in ( |B.6| )). It is obvious that 1 + K \s nonvanishing, 
so it can be reversed. Call A the inverse: 



1 



l + K{k) 



The inverse Fourier transform of A, indicated as A(x), is a distribution, because A(A;) 
k — > iboo. The following obvious equation holds: 



dx /\{x) {5{6 -x) + K{'Q - x)) = S{9) 



1 for 



and this allows to invert the convolution in ( |3.17D . Calling globally the right hand side of 

( p.l7| ) the inversion is: 



Z'^{^)= j dxA{^-x)J^{x). 



(3.18) 



The various terms in give different contributions, that shall be computed in the following. 
For the two terms ± 0, -) the convolution can be completely performed using ( A. 3 ): 



J dxA{^-x) (/.'(x±G,^; 



1 



2' cosh(i?±e) 

The effect of the inversion on (p'i'd, 1) gives an important object, that will be called G : 

4>'{x,i) 



dxAii) - x) 



2tt 



dxA{'d-x)K{x) = G{'&) 



(3.19) 



(3.20) 



Using Fourier transforms (see |A.3|) for A and (p' (B.6) the following expression can be obtained: 



GiO) 



p + 



27r 



+ lj 2tt ^ 'l + F 



dk e 



ike 



sinh 



tt{p — l)k 



K{k) 



(3.21) 



Tipk Trk 

2 sinh cosh — 

2 2 



The expression ( [B.6| ) used in the previous computation for the Fourier transform of (/)' holds only 
for 6 in the fundamental strip ([B.3D , i.e. for real solutions or close pairs, as in section |3.2| . This 
means that contributions to coming from wide roots require a different approach (see later). 
Instead the terms containing holes and close roots are completely arranged in this way. 
The Fourier transform of ( ^.21|) is obviously 



G{k) 



k{k) 
i + k{k)' 



This is an even function, real on the real axis {real analyticity) , positive if p > 1 and negative in the 
opposite case. The same properties, clearly, hold for G{9). The function G vanishes exponentially 
for large values of \k\: 



G 



~ e 
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(consequently also the function G{9) has a similar asymptotic behaviour). Because of this, the 
Fourier transform in (p. 21 ) at the points 



|9m0| = zbvr min(l,p) 

has a singularity. As before, this points correspond exactly to the limit between close roots and 
wide roots. This confirms exactly that wide roots require a different analysis. Call source the 



contribution to (|3.18|) by one wide root put in w: 



source = dx A(x) cft^-d — x — w,l) 



1 f f rlh 



X — w^l). 



Using the known expression (B^) for (j)' the computation can be performed into a closed form. The 
integral in x must be performed first. It is a Fourier transform, but because of the w contribution, 
it is different from the Fourier transform of 4>' calculated in ( [B.6D , that holds for ^ near the real 
axis. This is because ^mw is larger than the position of the first singularity of cp' , and the poles 
to take into account in the two cases are different. The result can be written in the following form: 



source = 2ttGii{'& — w) 
where the function Gn is defined by (for > 7rmin(l,p)): 



(3.22) 



27rGii{e) 
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coth — signQm(9) + coth 
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{iiT — 6 signSsm{9)) 



P 



sign9m(0) i 



1 



+ 



sinhS s'mh{9 — iTrpsignQm{9)) 



if p> 1 



if p < 1 



(3.23) 



A quite singular property, pointed out in ^j, is that this function can be interpreted as the so 
called second determination of the previous function G. The general definition of the second 
determination is: 



fiiiO) 



f{9)+f{e-iTTsignQm{e)) if p>l 
f{9) — f{9 — iTTpsignQm{9)) if p<l 



for 19^6*1 > 7rmin(l,p) (3.24) 



Applied to ( F21I) , it yields exactly ( ^2^) . 

Taking into account ( p. 19 , 3.20 , |3.22| ) in ( 3.18| ), the equation ( 3.17 ) takes the form of an 
integro-differential equation for Z^: 
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+- 

1 
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1 



cosh(i? + e) cosh(-i?-e) 

Mc Mw 

2vrG(i? - Cfc) - ^ 27:Gn{^ - Wk)+ 

k=l k=l 



Nh 

^InGi^-hk) 

k=l 



dxG{'d — X — ir])- 
dxG{'d — X + irj)- 
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^-lfe-iZN(x-iv)iz>^{x - irj) 

1 -|- (_l)<5g-i2jv(3;-«'7) 



(3.25) 
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This equation holds for both the cases of uj real or complex. The most common case is u; G M. 
This allows to write the last two lines in a more compact form, valid only if i} is on the real axis. 
The real analyticity of Zn and of G is used: 



J dxG{d 



iT])- 



1 -)- (^ — iyg,iZM(x+iri) 



if ^9 G 



Observe that the second factor in the integral terms seems to be the derivative of log [l + (— i)'^e^*'^^^^^*'^^] , 
but this substitution requires some care, because of the polydromy of the logarithmic function. 
Call, for the sake of simplicity, the argument of the log: 



1+ -ire 



5 iZNix+irj) 



(3.26) 



and use the fundamental branch as in the appendix |B|. If f{x) cross the cut of the log on the 
negative real axis, the function logpD f{x) has a jump by ±2Tri. The exact real point where 
this happens (suppose for simplicity only one such point; the extension is trivial) is given (see the 



figure 0)by the condition /(y^) G ] — cxd,0[, that is: 















V ' J 






(a) (b) 
Figure 3.3: Plane of the log. The circles have radius |e*^"(^+*'') I . In (a) it is < 1, in (b) it is > 1. 



1 + (_i)'5e'^^^^feT+^^) = and 9mZjv(yt + iv) < 0. (3.27) 
Only the expression 

logFD fix) - 2TTie{x - y|)e 
gives rise to a continuous function, e = ±1 is a sign that is positive if the cross is in the clockwise 



direction (from below to up) as in figure and negative in the opposite case. An important 
observation can be made by the following Taylor expansion: 



^eZ^iyi + irj) = ZN{y^) - —Z'lf{yi^) + ... 



(3.28) 



9m ZN{y^ +ir]) =7][ Z'j^{y^) - ^Z^'(y|) + ...). 



(3.29) 



It is clear that if and only if ry — > the condition ( p.27| ) becomes the definition of a special 
solution, i.e. y| is a special object only in the limit of r/ going to zero because in this case the 
second terms in ( 3.2j , 3.29D are negligible. Then the condition ( 3.27| ) becomes Z^(y|) < and 
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Z]si = 2ttI. The sign must be chosen positive: e = 1.^ For the general case r/ / 0, is a special 
object y shifted by a little bit. 

A completely analogous analysis can be performed for the lower integral, obtaining a y^. Prom 
(|3.28|) , because of the real analyticity of Zn, the following equation holds: = Vi = y- The result 
is: 

'-'''f:t:;:^^t"" -i-^^^ - (-d*.--'-") . - 0.30, 



Using this expression in ( 3.25 ), a new source term appears every times there is a such point ij: 

1 1 



+ 



Nh 

+ ^2^G(i?-/ifc) + 
fc=l 



cosh(T9 + e) cosh(i?-e) 

Mc Mw Ns 

2vrG(i? -Ck)-Y, '^^Giii'd - Wk) - 27r ^ (G(i? - ^ + iv) + G{d - yt - iv)) 

k=l k=l k=l 



(3.31) 



+- 



d 



5 ^iZ]^{x+iri) 



dxG{-d - X - ir])^\ogpj^ 1 + (-1) 
dx L 

dx L 



+ 



The sum on ijk has been indicated as sum on the specials, even if they are not exactly in the 
position of specials, because in the usual computations made on the NLIE rj takes a small value. 
The condition to be used for the correct position is ( 3.27 ). In the following, we omit the label 




Figure 3.4: Jump in the log: log/(y + v) — log/(y — ly) — 2-Ki 



FD, intending always the log in the fundamental branch. It is very important to note that special 
objects are not completely independent degrees of freedom, unlike the holes and complex solutions, 
that are fixed "a priori". Special objects appear when the derivative of Zat on a root or on a hole 
becomes negative. A special object then apparently appears two times in (|33|), both as a real 
root/hole and as a special root/hole. 



3.5 Non linear integral equation (II) 

From a Bethe Ansatz point of view, only the function Zj^ (and not the derivative) has a physical 
meaning, because it is required to obtain the position of various types of solution, given the 

In what follows, always this case will be assumed. 
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quantum numbers ( |3.2| ). To obtain Z^, an integration in 'd must be performed in ( p.31 ). All the 
source terms must be integrated. For the case of holes, specials and closes this is given by the 
function: 



xio) 



dx2'KG{x) 



(3.32) 



It is and odd and real analytic function. It asymptotic value is known by direct integration 

vr p — 1 



Xc 



x(+oo) 



2 p 



(3.33) 



Because G has always the same sign for every x, the x(^) is a monotonous function. 

The integration of the source for wide roots is simply the integral of ( |j.23 ). But this integration 
produces logarithms, then there is a problem of definition of the appropriate branch. The Zjv 
function is continuous in the analyticity strip, then the most correct choice for the wide sources 
satisfies this continuity. The forms that are proposed in the following do so (they have no jumps) 
for 9 with imaginary part 7rmax(l,p) > IQ'm^l > 7rmin(l,p) (the log is in the fundamental 
determination): 



xii{e) 



i sign(3'm0) ( logsinh ( - ) — logsinh 



ivr sign9m0~^ 



i sign($5m0) ( log ( — tanh - ) + log ( tanh ■ 



P 



inp sign9m0 



if p> 1 

if p < 1 

(3.34) 



The notation xii is used to indicate the integral of the Gjj. As in the previous case of G, this is 
the same as taking the second determination of x iii ( 3.32 ). The two forms corresponding to the 
opposite signs of the variable are related by the oddity: 

Xlli-G) = Xlli^) if \Qm0\ > 7rmin(l,p). 

This function can be continued to smallest values of the imaginary part, but obviously it has a 
line of discontinuity. The asymptotic limits are (now i? is real): 



lim Xll{'& ~ 

•d — > ±oo 



± (2Xoo - vr) if p > 1 

zb-TT if p <1 



(3.35) 



The final form obtained in this way is the so called fundamental non linear integral equation, also 
known as Destri-de Vega equation (see Q) for the counting function: 

smh "3 

Ziv(^9) = 2A^arctan + g{'d\'dj)+ 



cosh 



+ 



j ^Gi'd -x-irj) log (l + (_i)'5e^^iv(x+ir,)^ ^ 



(3.36) 



/ 



dx 



- I —Gi-d - X + ir/) log ( 1 + {-lye 



S „-iZN(x-ir]) 



+ a 
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where 

Nh Ns 

gi^l^k) = Y.xi'&- hk) - -yk + iv) + Xi^ -ilk- ir])) + 

k=l k=l 

(3.37) 

Mc Mw 
fc=l k=l 

Both the previous equations hold for 'd in the fundamental strip. Out of that, the analytic 
continuation of Zjv is required, because the first singularity of G is crossed. To obtain an equation 
expressing the analytic continuation of Z]y beyond this singularity, the second determination of 
the various terms appearing in the right hand side of ( p. 36 ) must be taken, as defined in (3.24). 
Notice that it is not the second determination of Zi\f but its analytic continuation: 



sinhi? 
arctan — — — ■ 
cosh B 



+ 



for |3=mi?| > min(l, p) : Zn{^) = 2N 
+g{^Wk)ii + I y G(t? -X- log (l + (-i)^e^^^(-+^^)) + (3.38) 

- j ^G{^-x + iv)iilog (l + (-l)^e-^^^("-^'')) +a// 

This explicit form of the analytic continuation requires to compute g{'&\'&k)ii- Remembering 
(|3.37|) , it is clear that for holes, specials and close roots the second determination of x{'& — '&j) (see 
( 3.34 )) appears. For the wide roots the second determination of xiii"^ — '^j) appears. It must be 
determined considering that 

Xw{'&,Wj) = xil{'& - Wj) 

is a function of 

X7/(i9 - Wj)ii = {xw{^, . (3.39) 

The explicit forms are not given there but they can be simply computed from the definition of 
second determination. 

The integral of the convolution term is possible because G vanishes exponentially at infinity. 
a is the integration constant and must be determined. 

There is a useful way to write the log term in ( |3.36| ), valid for real x and if r/ — > 0: 

Qn{x) = lim^^o 29m log (l + (_i)'5e»^iv(x+»»7)) = 

1 1 + tiYf,iZN{x+ir,) (3-40) 
= hm^^o - log ^ ^ ^_^y^.iz^^^_ir,) = (^^(^) + 27r 

The expression [A) mod 2tt means exactly — vr < A < tt. The second line in ( 3.40| ) holds only if 
Z'j^{x) > 0, because the branch cut of the log is not crossed and \Qm log/(x)| < vr (the notation 
■26|) is used). In the opposite case the branch cut is crossed and the previous expression can 



take values larger than vr: \Qm log/(x)| < 27r. In the limit case Z'j^{x) = (for example this 
happens when x — > oo) the values permitted are — vr < 9m log/(x) < vr. 
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At this point it is possible to explicitly calculate, from ( 3.36 ), the limit Zjv(+oo), using ( p. 36 
3^, |3]3|, |3^ ). The limit gives: 

Z^(+oo) = iVvr + xoo (Nh - 2Ns - Mc - 29{p - l)Mw) + 

(3.41 

+7rsign(p- 1)Mh/ + / dx G{x) Qn {+oo) + a. 



The limit QAr(+co) can be calculated by using in ( |O0| ) the limit of Zn (U). By a comparison 
of the ( 3.41D with the asymptotic values of Zn calculated in ( ^^ ) the value of a can be obtained: 



a 



p + l 

W hXoo 

p 



1 S 
2^ p+l 



+ 



I s 

2^ p+l 



UJ 

vr 



(3.42) 



In this expression the following term has been omitted: — 7rsign(p — 1) Msc- It is a multiple of vr 
and can be taken into account simply changing the value of b from to 1 (or viceversa) if Msc 
is odd, and shifting of an integer quantity the quantum numbers 

b={M + Msc)mod2 = {N-S + Msc)mod2 S {0, 1} . (3.43) 

This sort of manipulation on the quantum numbers has been explained at the end of the section 
Observe that the result is exactly a = if there is no twist. There is another important 



observation. In (|2.37| ) the twist term a; is invariant for the shift 

u; — > w + vr. 

The same sort of invariance is required in a for the NLIE (which is equivalent to Bethe equations) . 



It is simple to verify that the expression for a (|3.42D displays the following symmetry 

a — > a + 27r when u) lo + ti . 



Note that shifting a by 27r is an invariance of the NLIE ( 3.36| ), if an appropriate redefinition of 
the Bethe quantum numbers is made: Ij — > Ij + 1. This shift does not affects physical quantities, 
that depend only on the variables ^j. 

With the given value of the integration constant the equation ( 3.36| ) is complete. The quanti- 
zation condition (13.21) can now be written as: 



ZN{^j 



2-kI 



3 ' 



+ 



1 + 5 



(3.44) 



for the various solutions. For the wide roots remember the previous warning about the analytic 
continuation of the Z^ function. 

The previous equation is completely equivalent to Bethe equations. No new physics has been 
introduced, until this point. Simply an equation that generates the counting function Z^r has 
been obtained (the NLIE). 



3.6 Energy and momentum 



The expression (|2.4l|) can be written using the function 4> as: 

1 

^ = - E ("^(Q - ^i' V2) + 0(e + ^j, 1/2) - 27r) 



M 



(3.45) 



P=-Y1 - ^i' 1/2) - 4>{Q + ^j, 1/2)) + 2u; 
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The choice of the logarithmic branch in the energy ensures that the contribution of each real root 
is negative definite. This is consistent with the known ground state structure (as in section ^.5D , 
that is given by the maximal number of real roots. It will be clear that excitations give only 
positive contributions. 

It is possible to relate the previous expressions to the counting function. To do that, consider 
first the following quantity: 



M 



(3.46) 



It can be integrated to obtain the pieces appearing in ( ^.45 ). Clearly the following expressions 
hold: 



^0(9 -i?j,l/2) = / dxW{x 

j^l J asymp 



M M 

+ ^j, 1/2) = - E - V2) 

j=i j=i 



(3.47) 



dx W{x) 



asymp 



The integration requires to fix a constant. This can be simply done by computing for Q +oo 
the left hand side in (|3.47|) and imposing the equality with the primitive of W at the same limit. 
This has been indicated by the symbol 



asymp 



The function W admits an expression in terms of the counting function. The sum over roots in 
(|3.46|) can be expressed as in ( |3.8|) and using the same notations: 



Mr+Nh 



J' 



Nh Mc+Mw 
l/2)-Ycp'{e-hj,l/2)+ J2 </''(^-Ci,l/2). 



Using now the same trick used for Z'j^, as in ( 3.13 , 3.14 , 3.15|) , and the expression (|3.30| ) generating 
the "special" contribution, the following expression yields: 



Wi^) = I ^/{^ -x + i7j, 1/2)^ logpj, [l + (_i)5e-^^(--^) 



+ 



- X - ir/, 1/2)^ log^^ [l + (_i)'5e*^iv{x+..) 



+ 



Mc 



- J2 ^'{^ - hk, 1/2) + Y ^'(^ - Cfc' 1/2) + Y ^'(^ - ^fc' V2) + 

k=l k=l k=l 

Y -yk + ^V, 1/2) + 4>'{^ -Vk- iv, 1/2)) + J ^^'{^ - X, 1/2) Z'^ix). 



(3.48) 



Now the equation (|3.31|) can be put in the last integral in the previous expression. The terms 
"similar" can be collected (i.e. holes with holes, close roots with close roots, and so on) to obtain 
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the following form: 



N„ 

E 

fc=i 



(/>'(i?- /ifc,l/2) - / dx (p' -X, 1/2) G{x-hk) 



+ 



Mc 



k=l 

Mw 

+E 



<P'ii}-Ck,l/2)- / dx(P'{^ -X, 1/2) Gix-Ck) 



+ 



k=i 



(P'{^-Wk,l/2)- j dxcl)'{^ -X, 1/2) Gii{x-Wk) 

Ns 

+ J2 ["^'(^ -y^ + V2) + (/-'(i? - yfc - if], 1/2) + 

k=l 

dx cp'{^ - X, 1/2) -yk + iri) + G{x - yk - ir])) 

1 1 " 



+ 



+ 



29m / ^</>'(i9 - a; - ir/, 1/2) / dy [5{x - y) - G{x - y)] log^^ 1 + {-if e'^^^y^''^^ 



+ j ^ct>'{d-x,l/2)N 



27r 



cosh(i9 + 0) ~^ cosh('!? - 9) 
d 



+ 



dy 

Wh + Wc + Ww + Ws + Wu^ik + Wi 



where the notation introduced on the last line is obvious (the meaning of bulk will be clear in the 
following; label I is for interaction term or, that is the same, for integral). 
As in section 3^, there are different terms to analyze. 

First, rearranging a little bit the terms and introducing a delta-function, it is possible to see 
that in the case of holes, closes, specials and in the integral term there is a contribution by: 

-y) - G{x -y) = A(j; - y) 

(this equality is a simple consequence of the definition of A). The corresponding terms contain 
the integral: 



dx (j)'{x -^j,l/2)A{^ - x) 



1 



cosh(i9 — -dj) 

that is well known because it has been used in ( 3.19| ). This holds in the case of holes and closes 
and i?j is the corresponding rapidity. For the integral term and the specials the contribution is 
exactly: 

1 

cosh(i? — y — irj) 

where y is the integration variable or the special rapidity. 

The wide computation is more involved. The expression to be computed is 

Mw r „ 

WwW = Y. ^'{^-WkA/'i)- / dxc^'i^ -x,l/2) Gn{x-Wk) 
k=i^ 

and the most convenient way to do that is to use the Fourier transformation of all the terms. For 
— Wk, 1/2) and Gjj it has been used in the computation for ( p.23| ). For (j)' on the real axis it 
is in (|B.6D . The result is 



Wwi^) 



1 



cosh('!?) 



for p > 1 
+ 0(i?) for p < 1 
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In the case of p > 1 it simply cancels out. In the other case, there is no such cancellation, for 
generic wide positions. The explicit form of the contribution 0(t?) is quite long and shall not be 
written out at this point. The important fact is that it is a fast vanishing contribution for large 



■d. Observe that the second determination (|3.24|) of a whatever trigonometric hyperbolic function, 
for p > 1, is exactly zero, then the notation [1/ cosh(i?)]7/ will be used also in this case. 
The final form is: 



N 



H 



1 



cosh(i? - hk) 

Mc 



+ 



Ns 

E 

k=l 



1 



1 



cosh{'d - Ck) 
dx 

— - X, l/2)iV 



cosh('!? 



k=l 



Vk 

1 



if]) 



—2Qm 



dx 



1 



cosh(T? 

1 



Wk 



+ 



+ 



1 



cosh(a; + 0) cosh(3; — 0) 



cosh(i? — ijk + if]) 
+ 0{d)+ 

+ 



+ 



d 



log 



FD 



l + (-l)''e^ 



5 iZM{x+iri) 



2iT cosh('i? — X — irj) dx 

Integrating this function, as in (|3.47|) , the energy and momentum of this lattice system can be 
obtained. In all the terms appears 1/ coshx; its primitive is 



dx 



arctansinhx = 2 arctan tanh(x/2). 



coshx 

The last form is the most convenient, to calculate the continuum limit. Then: 



(3.49) 



E±P 



Nh 



Mc 



a- 



2 arctan tanh((e =F hk) /2) + ^ 2 arctan tanh((e =F Ck) /2) + 



Ns 



k=l 



k=l 



+ 2 (arctan tanh((0 ^ jjk — ^??)/2) + arctan tanh((9 =F yfc + ^^)/2)) + 



k=l 

Mw ^ 
+ ^ 2 arctan tanh((e =F Wk) /2)// + 0(6) =F / 
k=i 

f dx r 1 

+ / — <A(eTx,i/2)7v , . . + 



dx 



(3.50) 



2ir cosh(ib0 — X 
1 



■QAr(x) + S'7r+ 



cosh(x + 0) cosh(x — 0) 



±u- Ntt 



The term 0(0) vanishes for large 0. 

The function ( ^lo|) has been used, by putting — > 0, and observing that x is real. 



3.7 Continuum limit 

At the conclusion of the previous long computations there is an important observation: in the 
final expressions for energy and momentum there is no explicit contribution from the real roots. 
This also happens in the NLIE ( ^.36 ). On the lattice it is not a relevant observation, but in this 
section it shall be shown that it can justify a continuum limit procedure, — > oo and a — > 0, 
and a particle interpretation. In the lattice energy expression there is a bulk term (see later) i.e. 
a term increasing with N and completely independent from the solutions 'dj of Bethe equations. 
This bulk term can be subtracted and the remaining terms (due only to holes and complex roots) 
can be chosen in such a way that they describe particle excitations. 

To do that, in the limit procedure, they must be chosen in finite number (order 0(1)) and 
also its "rapidity" 'dj must be finite. Prom (3^) this means that only real roots can appear in the 
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asymptotic tails of Zjy. All the other lattice states must be discarded (as usual, lattice theory 
contains more states than the continuum theory). 

A state where holes and complex roots are not considered, is the hamiltonian vacuum. If they 
are considered, they give positive contribution to the energy, behaving as particle excitations on 
a vacuum state. Real roots are completely disappeared, in this limit. They can be interpreted 
(in the limit procedure) as a sort of Dirac sea that is the hamiltonian vacuum on with holes and 
complex roots built particle excitations. 

This is the structure required to have a consistent quantum field theory. 

All this interpretation is possible only because of the antiferromagnetic vacuum choice, made 
in section (p.5p. 

What will be proved now is that there is a consistent way to do continuum limit in the NLIE 
and in the expression of the energy. 



In it was shown that the correct way to do the continuum limit for this model on the 



light-cone is to send A^, G ^ oo connected in the following way: 



4N 

e»i„g— ^ (3.51) 



L = N a is the spatial dimension of the lattice (as in section |2.lD and stay fixed in the limit. M. 



is the renormalized physical mass. As in section |2.5| , the light-cone lattice is periodic in space 
direction. Then its continuum limit yields a compactified space of length L (space-time is a 
cylinder of circumference L). The lattice 6 vertex model becomes a field theory defined on this 
cylinder. The interpretation of what field theory is defined by this continuum limit, at the various 
values of the twist, is the job of the next chapter. 

To do the limit, consider first the counting equation. It takes the following form, on the 
continuum 

Nh - 2Ns = 2S + Mc + 2e{p- 1) Mw (3.52) 

because the structures in the tails are extremely simplified. The other terms appearing in the 
lattice counting equation (3J) in fact came from the configurations where in the tails there are 



special root/holes or self-conjugate of the so called first class (see the original paper |g||). Such 
configurations completely disappear in the continuum limit. 



The principal interest is to obtain the limit of the energy and momentum eigenvalues (|3.50D . 
Clearly, to fix the positions of the roots 1?^ the limit of ( 3.44 ) is required. This can be made using 
a continuum counting function that is defined by: 

Z{^) = lim Zn{'&). (3.53) 

A'' ^ oo 

This limit computation is made on the sequence of counting functions that are implicitly defined 
m (|]36|) at the various values of A^; the (|3.5l|) must be also taken onto account. The interesting 



fact is that this limit can be done explicitly in all the terms appearing in ( ^.36|) , i.e. a continuum 
NLIE can be obtained. To show that, the various terms are analyzed. 

The first term and the integral in ( 3.36| ) are trivial matter to compute. The source terms 



for holes, complex solutions and specials don't have any explicit dependence from A^. Only the 
positions that appear in this terms must be determined, obviously, by the continuum NLIE. Then 
the term g{'i}\i!}j) is unchanged. The result is: 



Z {-d) = ML sinh^ + g{'d\^j)+ 
+ I ^G{^ -x-irj) log (l + (-l)^e*^(^+^'')) + 

^G{'& -X + iri) log (l + ^-l)S^-^Z{x-^r^)\J ^ ^ 



(3.54) 
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The quantization conditions ( p.2| , 3.4^ ) can now be written as: 

1 + 5 



G Z + 



(3.55) 



for the various solutions. This terminates the hmit procedure on NLIE. Observe that on the con- 
tinuum only the implicit "definition" of Z by ( 3.54 ) is available, instead of and explicit expressions 
as in the case ( |3.1j ). 

Consider now in the energy expression (|3.50| ) the terms with the form 2/aarctan tanh((0 — 
9)/ 2). Using the following hyperbolic trigonometric identity 

vr _ 
arctan tanh(x) = — — arctan e ^ 

the following asymptotic behaviour can be obtained: 

^ 2 arctan tanh(G - 6*) « ^ - ]-M 

Then, collecting in the expression for [E it P)/2 the contribution of this type, yields (remember 
that there is a 57r in (IS^SOp): 



vr 



(2S -Nh + 2Ns + Mc + 2e{p - l)Mw) = 



(3.56) 



(the ( ^.52 ) has been used). 

One of the two integral contributions is 



nm =F — 

N~*QO L 



dx 



2tt cosh(ib0 — x) 



Qn{x) 



M f dx 
2^' 



„±x 



Q{x). 



The final result has been obtained by exchanging the limit and the integral. 

The other one is the last line in ( 3.50|) . It can be handled by shifting the integration variable 
and observing that one term in the square brackets is odd and gives a vanishing integral: 



^N, bulk 



iV2 

L 



I ^HQtx,1/2) 

/dx 
2^ 



N^ 
L 



1 



cosh(x + O) 
6(x,l/2) 



+ 



cosh(2e 



cosh(x — Q) 



vr 



IT 



This energy contribution is "source independent". It has been labeled by "bulk" because it will be 
shown that it diverges as N. It is interesting to compute the contribution to E and P, instead of 
the diverging contribution to {E ± P)/2: 



E 



Nbulk 



^N, hulk + ^N, hulk 



dx 



cosh(20 — a 



27r 



Ar2 

- — (7r + 7) 



P, 



Nhulk 



bulk 



E 



Nbulk 



2uj 



The first integral can be computed observing that (/> is a constant for very large values of x, and can 
be put out of the integral sign. The second one contains the difference of two equal contributions 
and only lo is not deleted. 
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Collecting all the terms follows: 

'Nh Ns 



E±P M 



(3.57) 



2 2 , 

\fc=i fc=i 

Mc Mw 



k=l k=l 
or, subtracting all the diverging contributions, 



± 

Nbulk 



/Nh Ns 

E = M ["Yj ~ E i^'^^^iVk + ^"U) + cosh(yfc -«??)) + 

\fc=l fc=l /o c;o\ 

Mc Mw ^ \ '-"^•^^^ 

— ^ cosh Cfc — ^ cosh// ~ 2~ sinh xQ(x) j 



k=l k=l 



and 



P = Y^ sinh /ifc — Y^ (sinh(yfc + irj) + sinh(y/s — irj)) + 



mT' Mw'=' ^ - N (3.59) 



sinh Cfc — sinh// ''^k + J cosh 2;Q(x) 



k=l k=l 



Observe that for p > 1 the second determination of hyperbolic functions vanishes. This means 
that wide roots contribute to the energy only in an implicit way, because they contribute to the 
position of other objects, by (|3.54| ). 

3.8 Physical interpretation 

The limit procedure described in the previous section is mathematically consistent, but the ques- 
tion is if from the physical point of view it describes a consistent quantum theory and allows for 
a meaningful physical interpretation. 

Before to propose it, an important remark must be made about the allowed values for 5. It is 
clear from ( 2.40| ) that on the lattice only integer and nonnegative values can be taken into account 



for S. But on the continuum the definition of S is no more related to the Bethe state (that is 



undefined), instead it is given implicitly by (3.52). Then, "a priori", there are no arguments that 



fixes its values to be integers. As our group showed in [g^], the half-integer choice for S is necessary 
(and gives completely consistent results) to describe odd numbers of particles. 

At this point the following physical interpretation can be proposed. It will be refined to 
describe the correspondence with particles. Also it will be supported by many arguments that 
will be clarified in the next chapters. 

Physical interpretation: 



the physical vacuum (hamiltonian ground state) corresponds to absence of sources (i.e. holes, 
complex); all the sources are excitations on this vacuum 

for u = and at the various values of 5 this theory describes the sine-Gordon/massive 
Thirring model on a finite space of size L; 2S is the topological charge and can take non- 
negative integer values. 
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• for the values 

0; = ^, k = l,...,q'-l (3.60) 

it describes the quantum reduction of sine-Gordon model, i.e. the massive integrable theory 
obtained perturbing the minimal model Vir{r,s) with the operator ^(1,3). 

Observe that it has been assumed that only nonnegative values of 5" are required to describe the 
whole Hilbert space of the theory. Indeed the theory is assumed charge-conjugation invariant then 
negative topological charge states have the same energy and momentum as their charge conjugate 
states. The assumption that all the states can be described by the NLIE is absolutely not trivial, 
or better still until this moment it is not available a general proof of this fact, but only a number 
of specific cases supports this conjecture. 

All this things are the argument of the next chapters. 



Chapter 4 

ANALYSIS OF THE CONTINUUM 
THEORY 

4.1 Principal questions 

In the next sections, a carefull analysis of the vacuum and of some excited states from NLIE will be 
performed at the various scales (various L), to completely understand the physical interpretation 
suggested at the end of the section |3]^. The principal question is what theory is described by 
equations ( 3.54 , 3.58 , 3.59| ). 



As scale parameter, as appears in ( |3.54| ), can be chosen equivalently the size L or the adimen- 
sional size I = where M. plays the role of a mass scale. The limit of very large / can be 

interpreted both as large size (that reproduces infinite Minkowski space-time) or large mass scale, 
that is an infrared point in a renormalization flow (IR). Finite size effects do not appear, in that 
limit. 

At the opposite limit of small I, important finite size effects are expected; a small value of I 
can be obtained with a small mass A^, then this case is an ultraviolet point in a renormalization 
flow (UV). The complete control on the scaling functions can be obtained analyzing the whole 
range of values / > 0. 



4.2 Connection with sine- Gordon/massive Thirring 



The section |3J| ends with a conjecture about the physical meaning of the model described by 
the continuum NLIE and the corresponding energy and momentum expressions. There are many 
arguments to support this interpretation. Some of them came from the well known properties of 
the 6- vertex Bethe Ansatz, other from the analysis of the NLIE itself. 

Consider the R matrix in ( 2.33 , 2.35|) . It has the same entries as the Rsg matrix that appears 



by putting sine-Gordon on a lattice, and obtaining the corresponding Bethe Ansatz equations 
(lattice Thirring B.A. gives exactly the same). The equality of R matrix and of Bethe Ansatz 
equations means that there is a common integrable structure in the two systems. 

Moreover, the lattice sine-Gordon B.A. can be obtained as the continuum limit — > oo (and 

4N r— , 

with 6 = log-——) of the inhomogeneous 6- vertex Bethe equations ( p.37D . The Bethe Ansatz 

Jvllj 

techniques, then, strongly suggest the previous interpretation. 

The 6- vertex model in its thermodynamics limit is critical, as shown in Q, then there is a 
conformal field theory describing this critical point. Prom the eigenvalues of the transfer matrix 
(obtained by B.A.) it is possible to extract the conformal properties, as shown in p3]. As shown 
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in |£3|, for 6- vertex model they reproduce exactly the equation ( |1.4| ) for u; = (that is the UV 
structure of sine-Gordon and massive Thirring) and the minimal models conformal weights for a; 
chosen as in ( p. 60 ). 

It has been shown, in that a fermion satisfying massive Thirring equations of motion can be 



constructed, as described in section (2A). It describes the scaling behaviour of the inhomogeneous 
6- vertex light-cone lattice dynamics. A similar argument can be made on the XXZ chain quantum 



hamiltonian, whose exponential is the transfer matrix of 6- vertex model. In |10| a (classical) field 
satisfying sine-Gordon equations is built from this XXZ chain. 

There is one more suggestion, that comes from the NLIE itself: the x function obtained in 
the derivation of the NLIE (3.32) is well known in literature, because it is exactly the logarithm 
of the soliton-soliton sine-Gordon S matrix: 



x(^?) = -ilogS++(^) 



(4.1) 



as m 



This arguments are the starting points for the analysis of NLIE itself that will be performed 
in the next sections. They completely confirm the physical interpretation given in section |3.7|. 



4.3 Some general facts about the IR limit 

This IR analysis is important, because it can connect NLIE with the (minkowskian) scattering 
theory. The limit L — > oo is not at all a thermodynamic limit, because the size of the system is 
scaled, but the number of particles is fixed. Then, the interaction among the physical particles is 
expected to cease affecting the energy and momentum, because their density is vanishing. Hence, 
E and P should approach finite limits equal to a free massive spectrum. In fact, for large I, the 
dominant term in ( 3.54| ) is the / sinhi?. Using it as first iteration in the convolution terms of ( ^.54 



I5^,g^), they become small (of order O (e ')) and can be dropped. The surviving part in NLIE 



can be seen as a dressed Bethe Ansatz giving constraints on the asymptotic states: 

Z{i}) = lsmhi} + g{i}\i}k) + a for / ^ oo. (4.2) 

Out of the fundamental strip where this equation holds, a similar expression can be written, using 
the second determination. Assume for the moment a = 0. Taking the exponential of the previous 
equation and imposing the quantization condition yields: 

gi«sinhi?jgig(i3j|i?fc) _ _|_]_ 

that has a structure similar to the usual quantization condition of particles in a box: 



e-^--'l[S{Xj-Xk) = ±l (4.3) 



{+ is referred to periodic and — to antiperiodic boundary conditions). The equation (O) is not a 
Bethe Ansatz equation]^ , in general, because the that appear in it are real rapidities of particles, 
instead of the rapidities of the pseudoparticles (??j) required to built Bethe states. But in (4.2) 



can appear complex "rapidities". The correct correspondence between the box quantization ( [4.3] ) 
and the infrared NLIE ( ^ ) will be carefully analyzed later. The important fact, for the moment, 
is that the comparison between this two equations allows to read out scattering data from NLIE 
and this can be compared with the known S-matrix of the models to which the NLIE is referred. 

^ It is called Dressed Bethe Ansatz because it contains real particles instead of pseudoparticles appearing in 
usual B.A. 
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Another important IR observation is that the derivative of the counting function is a very 
large positive number and no one special root/hole can takes place 

= / cosh ?? + £/' (i? I for I ^ oo. 

This suggest that special roots/holes are not physical particles. Instead, as it is clearly shown in 



section 3.4, they are a mathematical artifact of the logarithmic contribution that appears in the 



convolution term. 

The final important observation in IR analysis comes from quantization conditions ( 3.44 ) 

2TrIj = I sinh-iJj + g{-&j\^k) + a 

The number Ij is real. Distinguishing the real and imaginary part of this equation gives: 

27r/j = /sinh ^e'&j cos Qrwdj + SRe {g{'&j\'&k) + a) 
= / cosh 3f?e T?j sin Qm i?j + {g{'&j\'&k) + 

To deal with the first equation, the real part of the i)j for increasing values of the large scale I 
moves toward the origin as 



I cos 9m 

Prom the second equation, the imaginary part must develops a singularity from Qm {g{'i)j\i}k)) 
to annihilate the rapidly diverging term I sin ^rwdj. Because the positions of singularities are 
known, this simple argument can be used to fix (at the IR) the imaginary part of the sources. 
Being the real part zero, the source position at IR is completely fixed. 

4.4 The intermediate regions 

The NLIE (|3.54|) contains a source term g{'d\'dj) which is specified by giving the root/hole structure 
of the given state and depends on the positions of the holes, special roots/holes and complex 
roots. These positions in turn are fixed by the Bethe quantization conditions ( ^.55|) . The NLIE 
supplemented by the quantization conditions gives a set of coupled nonlinear equations in the 
function Z{'&) and the variables that can be solved numerically by an iterative procedure. First 
one chooses a starting position for the sources i^j. Then one iterates the integral equation (using 
fast Fourier transform to evaluate the convolution) to update the counting function Z. Using this 
new Z, an improved determination of position of the can be obtained, and is fed back into the 
integral equation for a new iteration cycle. The process is repeated until the solution is found 



to a prescribed precision (usually 10~^). As explained in section O, for large I the source term 
dominates, while the correction coming from the integral term is exponentially small. Therefore it 
is reasonable to expect that the further one goes to the IR regime the faster the iteration converges 
which is in fact what happens in the computations. Hence it is preferable to start iterating at the 
largest desired value of / and decrease the volume gradually, always taking as a starting point at 
the next value of the volume the solution found at the previous value. 

As a preliminar analysis, the general behaviour of the Z{x) function is that of sinhx for large 
I. Decreasing the scale, the tangent in the central portion of Z{x) for x G M becomes more 
and more horizontal; instead for large \x\ the leading term is always the hyperbolic sine. This 
central position broadens to a single or double plateau system that extends to all |a;|'s smaller than 
log(2//), and rapidly disappears for greater values, in favor of the dominant exponential growth. 
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It can happens, when decreasing Z, that special holes appear. But there is the so called number 
of effective holes, that is given by 

NH,eff = NH-2Ns (4.4) 
that is a constant independent of I. This is a consequence of one simple fact: when a globally 



increasing function (as Z is) makes a small oscillation as in figure 4J_ there are an odd number 
of points that intersect the horizontal line corresponding to a certain quantum number. Then to 




Figure 4.1: The typical behaviour of the Z function when zi is a special root /hole 

every real solution (when Zis increasing) correspond three (of five or seven or ...) real solutions 
in the case of specials. Of this odd number of solutions, only one is a real root (no Bethe roots 
can have the same quantum numbers). All the other must be holes (specials or normal). Then 



(as in figure 4.1) 



1 root for large / = 1 root + 2 holes — 2x1 special r. or h. for small I. 



The general case is ( |4.4| ). The physical meaning of the holes, then, must be related to the effective 
number N^^eff-, not to A'^f^. various 

For the finite values of the scale, it is interesting to made the comparison of the numerical data 
from the NLIE predictions with those obtained with the Truncated Conformal Space Approach 



(TCSA, see section ( |1.4D ) at several values of the parameter p. For illustration, some interesting 
cases will be presented. For the various values of the coupling constant, our group obtained a 
spectacular agreement between the results obtained by the two methods, up to deviations of order 
10"^ - 10-3 (see I, 1^, m). The deviation grows with the volume L, exactly as expected for 
truncation errors. As also told in (p^), in the attractive regime {p < 1) the TCSA convergence is 
fast. The opposite happens for the NLIE, i.e. it converges faster in the repulsive regime {p > 1) 
and for large volume L, but in general it is possible to have precisions of order 10"'' (and higher) 
in both the regimes. 

By studying other parameters such as the mass gap, the breather-soliton mass ratios and the 
rate of convergence of the energy levels with increasing the value of Ecut, the small differences 
between the two methods can be clearly attributed to the inaccuracy in the TCS data. 



4.5 The UV limit computation 



The ultraviolet limit of the states described by the NLIE will be examined, in order to compare 
it with the known facts about the UV limit of sG/mTh theory and conformal minimal models, 
outlined in sections (L2, 1.3). 
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More explicitly, it is known, from conformal perturbation theory, that the behaviour of energy 
and momentum for very small L are given by: 



E{L) = = (c - 12(A + A)) + . . . 

6L 6L ^ » (4.5) 

P(L) = -(A-A) + ... 

where the ellipsis are most regular terms in the space size L. Therefore, to determine the UV 

behaviour, in E and P only terms containing — must be retained. Consider, e.g., the hole 

contribution to the energy: M. cosh/i. In order to have this hole, this term must behaves as 1/L, 
that is possible only if 

2 

= finite ± logy. (4.6) 

This means that only holes with this rapidity can contribute to UV. The same argument applies 

to specials and complex solutions. 

The behaviour of the sources for / ^ can be classified by three possibilities: they can 

2 

remain finite (they are called central), or they can move towards the two infinities as itlogy 

[right/left-movers). The finite parts of their positions can be obtained by extracting the divergent 
part: 

= |^±±logy,^?o|. (4.7) 

(the I dependence of the roots has been made manifest). The number of right^left moving resp. 
central holes is indicated by N^' and similarly the numbers Ng ' , M^j' and M^/ are introduced. 



The finite parts in (|4.7D satisfy a modified version of the NLIE, known as kink equation. To obtain 
that, observe that Z has an implicit dependence from / that can be made manifest writing Z{X, I). 
Define the kink functions with the following expression: 

Z±(i?) = lini Z (^t?±logy,/^ (4.8) 

The term "kink" has its origin in the fact that the asymptotic form of the function Z{'i}) has 
two plateaus stretching between the central region and the regions of the left/right movers. The 
functions Z±{'d) describe the interpolation between the plateaus and the asymptotic behaviour of 
Z{'d) at the corresponding infinity. Note that if there are no central objects then the two plateaus 
merge into a single one stretching from the left movers' region to the right movers' one. Using 
( [4.7| , [4.8D in the NLIE ( 3.54| ), the source terms may behave in different ways. Consider e.g. one 



hole h{l) of the type /i^. Its source term for the kink "+" is: 

lini X ^-i? + log y - h{l)J = x('!9 - h'^). 
There are terms of this of type. For a h~ hole the source contribution to the same "+" kink is 
lim X fi? + logy - h{l)j = lim X f A - /i" + 21ogy j = Xoo 
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(Xoo is defined in ( 3.33 )). One hole behaves in the same way. Then there are Njj — = 
+ Njj such terms. Analogous arguments apply to "-" kink and to the other roots. For wide 

roots the limits of second determination are given in ( p.35| ). 
The following definitions are introduced, according to ( t3.52|) : 

5±'0 = i[iv±'0 - 2iV|'° - M±'° - 2M±'°0(p - 1)] 

Observe that S is always integer, whereas 5^'" can be half integer; moreover, in some cases, 
can be negative. All this numbers are interpreted as the spin of the right, left movers and fixed 
solutions. Clearly 

S = S+ + S' + 5°. 

Using all this computations in the definition (^!^), the NLIE for the "+" and "-" kinks reads: 



Z±(i?) = ±e^^ + a + g±{^)+ J dxG{'d - x)Q±{x) (4.9) 
where the following definitions are used: 



1? 



k=i 



k=l k=l k=l 

1% = ±sign(p - 1) ^ {Mw - M± ) 



(4.10) 



and the function Q±{x) is related to Z± as Z to Q in ( p. 40 ) 



This equations allow to write the quantization conditions as: 

Z± Uf] = 27r/±. (4.11) 



It is a matter of convenience to put the apex it on the quantum numbers. Indeed, they do not 
change in the limit procedure: they are exactly the same as for finite /. 

By simply taking the limit / — > in NLIE reads an equation for the fixed objects ??": 

Zo(i?) = lim Z(t?,/) = a + 5o(i?) + j dxG{{) - x)Qq{x) (4.12) 



where 



^ fe=i 

N°s M° Ml 

E (x(^ - fk + ^v) + -fk- iv)) - E - - E ^^^(^ - ; 

fc=l k=l k=l 

«o, = sign(p-l)i(M^-M+). 
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As in the previous case the function Qo{x) is related to the corresponding Zq{x) by the usual 
expression ( 3.40| ). 

In the following the asymptotic values of Z±^q and the corresponding values for Q will play an 
important role. Equations ( [4.9| ) simply gives 

Z+(+oo) = +00 = -Z_(-oo), (4.13) 

consequently the value for Q{x) can be obtained using ( |3.40| ): 



Q±(±cx)) = (4.14) 

(to compute it, before do the limit in x and later the limit in 77, because the NLIE holds with a 
finite rj). Moreover it can be shown that 

g+{-oo) = 5o(+oo) = 2xoo{S - 2S+) + 27r/c+ 
g-{+oo) = 5'o(-oo) = -2x00(5* - 2S') + 2iTk^ 



and 



"2 



= ±- sign {p - 1) {Mw - 2M^) = + integer. 



A direct consequence of this is that 

Z+(-cx)) = Zo(+oo) and Z_(+oo) = Zo(-oo) (4.15) 

because they satisfy the same equation. Remember that for small / the counting function gives, in 
the central region, one or two plateaus that stretch to the infinities. Then the previous equations 
simply mean that Zq{+oo) is the right plateau height and and Zo(— 00) the left plateau height. If 
there are no central objects, they take the same value, yielding so a one plateau system. 

The explicit values can be computed, using the kink equations. The integral can be handled 
because the Q in the asymptotic limit takes a constant value and the kernel G rapidly converges 
to zero: 

Z+(-oo) = a + g+{-(X)) + — Q+(-oo) 

(4.16) 

Z_(+oo) =Q + 5_(+CX)) + ^Q_(+00). 

TT 



They are called plateau equations. From the definition of Q in ( 3.40 ), and making attention to 



the warning concerning the range of values allowed for it, it follows that 

Z±(=Foo) = Q±(=Foo) +1:5 + 2Tik± 

with an appropriate choice of the integers k± such that the condition — vr < Q+(— c>o) < tt holds 
(remember that can be half-integer). The solution of the plateau equation is then 

w± = Q±(Too) = ±27r^^(S - 2S^) + 2^-^ (--5 + 2k%- 2k±) . (4.17) 

p + 1 p+lVvr J 

Notice that for p > 1 there can be cases without solution. If instead it is there, there is no 

p 

ambiguity in the choice of fc-t, because its contribution is a multiple of 47r > 27r, that is 

p + 1 

larger than the range of values allowed for Q±. Instead, for p < 1 the solution always exists, 
but it can happens that two or more choices of k± satisfy the condition on the range of Q±. It 
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will be clear later, that this phenomenon is related to the fact that wide roots become excitations 
independent of the other one's, which is manifestly evident from the fact that they no longer 
contribute to the spin of the state S and that the asymptotic value of their source contribution 



Xii in the attractive regime is simply a vr (see 3.35 ). 



Before to continue the general computation, the kink equations for the wide roots will be 
explicitly written. They can be obtained by using the same procedure used for (f4.9D , applied to 



the correct expression for Z^'d) given in (3.38) 



(4.18) 



k=l k=l 



- E ^^^(^ - 4) - E ^^^(^ - "^k)!! ± - l)4Xoo(S - 5±) + 27rL± (sign9mi?)+ 
fc=l k=l 

+ J dx Gn{-d — x)Q{x) + an for IQ'mi?! > min(l, p) 

where the notation introduced for L\y is so complicated because it depends on the sign of the 
imaginary part of 

Mw-M^ 

Ltr{signQm-d) = - lim [x//('i? - "Wfe + 2;)// =F 4(Xoo - vr)6l(p - 1)] = vr • integer 

(4.19) 



X — >±oo 

fc=l 



The explicit form is quite complicated. Then this contribution must be computed case by case. 
Observe that the term ef j vanishes for p > 1 and a// vanishes for p < 1 . 

Now, the ultraviolet limit on the energy and momentum expressions will be done. To this end, 
substitute in the energy and momentum expressions ( 3.5S , 3.59| ) the ultraviolet behaviour of the 



roots (|4.7|) and retain only terms in — . The kinetic terms give contributions as 

7WcoshT?fc ~ ye^'^fc . (4.20) 

The integral terms must be calculated separately. Their form in the energy and momentum 
expressions is, respectively: 

-M [ ^smhxQ{x) = -M [ ?^[e" Q(x) - e"" Q(x)] , 
J 27r J ZTT 2 

-M [ ^coshxQix) = -M [ ^^[e"Q(x)+e-"Q(x)], 

J 2lT J ZTT 2 

where the contribution depending on is called "+" kink contribution, and the other term is the 
kink contribution. Prom the definition itself of Q± observe that, in the limit of very small I, it 
is possible to write for Q: 

Q ± log y, A ~ Q± (x) + q±{x, I) , (4.21) 



where the two functions q±{x,l) vanishes in the I — > limit. The expression ( [1.21 ) can be 



substituted in the integral form for the two kinks. The integration variable can be shifted by 
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X x + log(2//) then 

where the symbol = means that only the terms of order 1/L are retained. Similarly for the "-" 
kink term 

At this point, it is possible to express energy and momentum in a way dependent only on quantities 
which are finite in the UV limit. Using then ( 4.2C , f4.22| , [4. 231) in (|4.5| ), and rearranging the 



expression in terms of A, A (i.e. restoring the light-cone coordinates) follows: 




As a notation, A = A"*" and A = A . The NLIE can be used now to express the various terms 



appearing in the sums and the integrals. Consider, for example, the quantization condition ( 4.11 ) 



for the holes hf. Using the oddity of the x(x) function, it can be arranged in a more convenient 



way: 



j=l j=l j=l \k=l 



+ E ^^^f - 4) + E ^^^^^f - ^fc ) - (« ± 2xoo(5 - 5±) + 27r/ 



w) 



k=l k=l 



^ 2 9m y yG(/i± -x-i7]) log (l + (-)5e*^±(^+^'')) 

(the nonzero r] has been restored for later convenience). For the other objects, similar forms hold 
(remember that the yk are defined by ( ^.27[ )). All the equations so obtained must be substituted 
in the expression for A^. Now, observe that all the sums of x{^), Xi'/(^) and xii{'^)ii cancel 
completely for the oddity of the functions. Moreover in the integral term the following substitution 
(for all the type of sources) can be made 

2TiG{hf -x) = x'{x-hf). (4.24) 

The result is: 



±£± T29m j ^^y,^(x + ir/)log(l + (-)V^±(-+*^)). 



58 



CHAPTER 4. ANALYSIS OF THE CONTINUUM THEORy 



The following notation has been introduced: 



iV± 



















j=i i=i i=i i=i 

The new constant appearing in that equation takes into account for wide roots, and is given by 
(see also (|4.18D ): 

£± = (sign$>m«;±) = -7rsign(p - 1)M± 2(5 - + 



fc=l 



EE lim Xii{w^ -Wk + x)n^A{xoo-T^)0{p-l) 



TT • integer 



It must be computed case by case. 

The following function has been introduced: 



k=l k=l 



(4.25) 



- E - - E ^^^(^ - ) = ^± w - / ^^^(^ - ^)2±(^) - « 

fc=l A:=l 

Now, with the help of the computations shown in appendix the integral becomes: 
and the UV limit computation ends with a close form for the conformal dimensions: 

vr ' lovr^ p 

This means that the UV limit admits an exact computation of the spectrum. The identification 
of the UV states is one of the fundamental steps to understand the physical interpretation of the 
continuum model so far defined. 

An interesting phenomenon is that the conformal weights obtained depend only on very generic 
features of the source configuration such as the asymptotics of the left and right moving sources, the 
total spin and the number of self- conjugate roots. This means that if a certain source configuration 
is given, one can add new sources separately to the right and the left moving part in such a way 
that they are separately neutral (i.e. do not change S"*" and S~ , the total spin and the number 
of self-conjugate roots). In this way the primary weights do not change, but generally the term 

— 21^ — — is increased and so descendents of the initial state are created. An example: 
states which have 5 = = Msc and = are all descendants of the vacuum, however 
complicated their actual source configurations are. 

It is important also to drive attention to the well known fact (see the section 3^) that (at 
the moment) there is no proof that the scaling functions obtained by the method of NLIE span 
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the complete space of states. This is an extremely difficult problem due to the following two 
circumstances: (1) the dependence of the UV conformal weights from the parameters of the 
source configuration is rather complicated and (2) to obtain the allowed values of the complex 
roots one has to carefully examine the IR limit as suggested in section O. It will be also clear 
that the same state can be realized by different root configurations depending on the regime (e.g. 
the {ss)± states: scattering of soliton and antisoliton in even and odd wave functions; see later). 



4.6 Sine-Gordon and massive Thirring 



As suggested in section (3.7), the choice a; = a = is supposed to describe the sine-Gordon and 
massive Thirring models on a cylinder. In this section the corresponding NLIE will be analyzed, 
starting with the vacuum state. Obviously, the coupling p is expected to be the same that appears 
in the s-G and mTh lagrangians via the equation (|1.12| ). 

Starting from a lattice with 2A'^ sites one finds that the antiferromagnetic ground state which 
has spin 5 = 0, when written in terms of Bethe vectors depends on M = roots of the Bethe 
equations, all of which are real, as indicated in ( |2.5D . This ground state is expected to correspond 
to the vacuum of the field theory. However, there are two possibilities: one for even and the 
other for A^ odd, corresponding in the continuum to the choices 6 = or 6 = 1. As it will be 
shown in the sequel, only one of these states can be identified with the vacuum for a local field 
theory having a c = 1 UV limiting CFT. The UV dimensions of the vacuum are = (because 



the theory is assumed unitary). Choosing 6 = 0, the expression ( 4.26 ) gives a value consistent 
with the interpretation proposed, i.e. 

c = l, 

as obtained in the paper |§. At the IR, because no holes neither complex roots are considered, 
the energy and momentum ( 3.5£ , 3.59| ) vanishes, that is what is expected for the vacuum. The 



numerical iteration of the NLIE gives the scaling energy shown in figure |4.2| . In that figure, 
the value of the predicted bulk energy ( |1.21| ) has been subtracted from TCSA data in order to 
normalize them at the same way as the NLIE data are. The agreement of the two data is to order 
10"^ - 10-3. 



As a consequence of this vacuum analysis, the equations for the asymptotic behaviour (|4.17 ) 



and for the conformal dimensions ([4.26| ) take a simpler form (a = and c = 1). If the physical 
interpretation is consistent, the expected conformal dimensions must have the form (|l.4|) , that is 
a sum of powers of R^ = (as is (|1.12|) ). Then it is convenient to express A='= in powers of 
using also the expression (4.17) for Q^(=Foo). The final form is 



= + ± +N± = 

2 ^ (4.27) 



p + 1' ^ ' 16 p 2 



± -mR^ \ + A^± 



2 V ^ 2 

where the following identifications have been made: 



m = 2S 

u± , . \^^c nc±^ (4.28) 



n± = - fc± + T{S- 25±) 



and 



N± = ± (4 - 2/± - /± - /± ) T 5±2/± ± £± - 2 {S^Y T 25± (^^-k^ + k^). (4.29) 
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Figure 4.2: The first few energy levels (TCSA) in the vacuum (m — 0) sector at p — j (plotted with 
diamonds) for Ecut — 17.0 (the dimension of the space is 4141) and the NLIE prediction for the vacuum 
scaling function (shown with a dotted line). 



The similarity with (L4) is manifest, 
required. The first one requires that 



To match an exact correspondence two conditions are 



\n+\ 



\n- 



(4.30) 



The second condition is that one must be able to chose the sign of n in such a way that the left 
and right descendent numbers are integer. Consider the case where n_|_ = n_: then (using ( 4. 27] )) 

follows that both N± must be integers. Instead, in the case n+ = — n_ a combination of N± and 

n±m 

± must be mteger. 



Since the contribution comes exclusively from the last term in ([4.26| ), the one-plateau 
systems (i.e. when Q_|_(— oo) = Q_(+oo)) satisfy trivially the condition (|4.30D , but even in that 
case it is not clear whether the descendent numbers are integers. At the time of this writing 
the analysis of our group cannot exclude the possibility that there exist some states for which 
the conformal weights cannot be interpreted within the framework of c = 1 CFT. However, in 
the numerous explicit examples we have calculated so far we have not found any configuration of 
sources which lead to such behaviour. If it happens for some configuration of the sources then the 
corresponding scaling function must be excluded from the spectrum. 

From the second line of ( [4.28| ) one obtains a useful relation (m = 25" G Z) 



n± 



m + 6 + Msg 



+ integers. 



Using now the classification of UV conformal operators for s-G and mTh in the two algebras .4;, 
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and Af given in section 1.3, one can obtain the following rule: 

m + 6 + Msc ^ ^ . ■ r A 

— G & tor sine Gordon 



6 + M^c ... (4-31) 
G Z for massive Thirring. 



With reference to figure 1.2, from this rule follows that all the four sectors can be accessed by 
NLIE, because the sectors I and II describe sine-Gordon, I and III describe massive Thirring and 

S + Msc 

the sector IV can be obtained by m even and half-integer. If there are no self-conjugate 

roots, the rule simplifies: 5 = for mTh (I and III) and (5 = 1 for the sectors II and IV. The 
sector IV, that also is accessible by NLIE, contains non local states. Then NLIE describes also 
this non local sector. 

4.6.1 "-" vacuum 

Consider the case of g {'&{'& j) = (i.e. no sources). There are two possibilities 6 = 0, 1. One is the 
vacuum, considered previously, the other is (5 = 1. The (|4.17 ) yields 



Q±(Too) = -27r^(l + 2/c±), (4.32) 
p+l 

which admits solution only in the attractive regime when 

|Q+(-oo)| = |Q_(+oo)| =2^-^. (4.33) 

The value is not determined uniquely due to the contribution of wide roots, but it can be fixed 
using information from the repulsive regime. 



In that case, by performing the usual iteration procedure described in section [4.4 one finds 
that only for I > Iq the iteration converges. For smaller values of I there is no convergence. The 
"breakdown" volume Iq depends hy p: for p = 1.5 it is Iq = 0.3 (no breakdown in attractive regime). 
What happens is that the real root at the origin is a special one with Z'{0) < 0. By ( p.52| ), because 
S = 0, also two holes must appear. All of which sources have Bethe quantum number zero. One 

of the holes moves to the left and the other one to the right and so = -, while the special root 

is central (the configuration is completely symmetric, then Z{—x) = —Z{x)). This allows for a 
unique solution of the plateau equations: 

Q+{-oo) = ±—-, Z+(-oo) = -Z_(+oo)=7r^<0. 
p+l p+l 

Observe that Z_(+oo) > Z+(— oo). This compared with ( [I.ISP means that there are two plateaus 
and the left is higher than the right, i.e. Z'{0) < 0. The picture is completely consistent. This 
phenomenon can be thought of as the splitting of the root at the origin into a special root and 
two holes and is one of the general mechanisms in which the special sources are generated (see 



figure [1.31 ). Another mechanism will be explored when the UV limit of the soliton-antisoliton 
states will be examined. Let us suppose that the counting function Z is analytic as a function of p 
which is plausible because all terms in the NLIE for this state are analytic in the coupling. This 
determines which branch to choose for the plateau values in the attractive regime. One obtains 
the final result (for both the attractive and the repulsive regime) 

4p+l 8i?2 ^±1/2,0' 

which are the conformal weights of the vertex operators V(_|_i/2,o)- The actual UV limit can be a 
linear combination of these operators; in any case, it is not contained in the UV spectrum of the 
sG/mTh theory, but is a state in the nonlocal sector (sector IV in figure |T3 
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a) b) 

Figure 4.3: The UV behaviour of Z(two plateaus system) for (a) p < 1 and (b) p > 1. In (b) the root 
at the origin turned into a special one and two normal holes appeared at the points where Z crosses the 
horizontal axis with a positive derivative. 



4.6.2 Pure hole states 

The analysis of pure holes states is the simplest one. In the IR limit one such state gives (using 
(PI): 

Nh 

Z(i9) = /sinhi? + ^xi^- hj) , Z{hj) = 2^1 j. (4.34) 
j=i 

The holes "rapidities" are real. Using the observation ( ^j] ) the equations ( |4.34 ) show the same 
structure of the box quantization equations for physical particles (|4.3|) . Then the holes can be 
interpreted as particles with rapidity hj, i.e. they are the solitons of sine-Gordon. Pure holes 
states are states with only solitons. 

The IR limit puts no restriction on the quantum numbers Ik of the holes; instead the rapidities 
converge to zero. This seems to be a very general feature that remains true even when complex 
roots are allowed. 

The UV calculations, for one hole, give for 5 = 1 (the hole can be put fixed in the origin): 



S=-- 5± = 0. 
2' 

Solving the equation for Q±(=Fc«) gives: 

m = 1 

n-i- = n_ = 

that correspond to the operators Vb,±i. They create the sine-Gordon soliton (see also ||2^). For 
one hole with 5 = 0, using 



S+ = -; S- = S = {), 
2' 



the UV computation gives the dimensions of 1 or 1^-1/2, -1 that describe the Thirring fermion. 
The UV picture is consistent with the IR. The behaviour at intermediate values of I is shown in 



figure |4.4|where it appears that the interpretation of the hole with the soliton in consistent. This is 
the simplest case where the NLIE was used, in |28], for an odd number of particles. Table 4.1 gives 



4.6. SINE-GORDON AND MASSIVE THIRRING 



63 



^ ^ ^ ^---^---^ o ^ o ^ o <> 



01 23456789 10 

I 

Figure 4.4: Comparison of tiie results for the one-soliton state coming NLIE and TCSA at p — |. 
Continuous line is the NLIE result, while the diamonds represent the TCSA data. 
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Table 4.1: Numeric al co mparison of the energy levels predicted by the NLIE to the TCS data for the 
case depicted in figure iA. The NLIE data are exact to the precision shown. 
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an idea about the numerical magnitude of the difference. Note that the deviations are extremely 
small for small values of / and they grow with the volume, exactly as expected for truncation 
errors. 

For two solitons case the two possibilities, corresponding to (5 = 0, 1 must be checked. The 
simplest cases are: 

1. Two holes quantised with Ii^2 = S = 0. It yields: 

corresponding to the operator V(o,2) > which is the UV limit of the lowest-lying two-soliton 
state, as it can also be seen from the TCS data. If instead of the minimal choice Ii^2 we 
take a nonminimal one /+ = 3/2, 5/2, /„ = —3/2, —5/2, we obtain 

which corresponds to descendents of V(o,2)- This is a general phenomenon: the "minimal" 
choice of quantum numbers yields the primary state, while the nonminimal choices give rise 
to descendents. 

2. Two holes quantized with Ii^2 = il, 5 = 1, as proposed in ||8|. In the repulsive regime 
a special root is required, as in the case of the 6 = 1 vacuum state. The result (in both 
regimes) is: 

This state is a linear combination of a-iV(i/2,2) and a_iV(_i/2,2) which means that it is not 
contained in the local operator algebras of sG/mTh theories. 

3. Two holes quantized with Ii = 0, I2 = ±1, 6=1. Consider in detail the case of I2 = +1 
since the other one is similar. Suppose that the hole with /i = is a left mover and the 
other one is a right mover (the other possibilities lead to a contradiction). We find a solution 
to the plateau equation only in the attractive regime: 

Q+(-oo) = ±27r^-, Q_(+oo) = ±27r^- . 

p+1 p+1 

In the repulsive regime the hole with quantum number becomes a special hole y and emits 
other two ordinary holes each quantised with zero. We obtain for the plateau values 

Q+(-oo) = Q-(+oo) = 27r^— . 

p+1 

The conformal weights turn out to be: 

1 1 

A+ = IH ; = A+ „ „ , 

4p(p+l) V2,2' 

A- = - i = A+ 

4p{p+l) 1/2,2- 

These are the conformal weights of the vertex operator V(i/2,2)- Performing a similar calcu- 
lation for I2 = —1 we obtain the weights of V(„i/2,2)- 
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The calculation for a generic number of holes proceeds in complete analogy with the cases 
treated until now, we only give a summary of the results. The conformal families which we obtain 
depend on how many holes move to the left and to the right. The primaries are obtained for 
the minimal choice of the hole quantum numbers; increasing the quantum numbers we obtain 
secondary states, as it was pointed out in the case with two holes. The states we obtain are in 
the conformal family of a vertex operator Vn,m with 



m 



NH,eff , n G Z + 



As a consequence, all states with 6 = are contained in the UV spectrum of the sG/mTh theory 
while the ones with 6 = 1 are not, in agreement with the relation ( [4.31| ). The complete expression 
for n is somewhat complicated, but we give it in the case of symmetric {N^ = Njj = 
configurations: 



n = 



5 = 0, 



n 



The comparison with TCSA completely confirms the analysis made until now. The example 



of four solitons states in repulsive regime is in figure 4. 5. As usual in TCSA, one can see that 

50 




Figure 4.5: The first few energy levels (TCSA) for m = A solitons at p = 1.5 (plotted with diamonds) 
for Ecut = 22.5 (the dimension of the space is 4149) and the NLIE prediction for the four holes scaling 
functions with quantum numbers (— |, —\, 5, |), (— |, —5, 5, |)and (— |, — |, |, |) (shown wit continuous 
lines). 



the truncation errors become larger; at values of the volume I close to 10 the deviation can be 
observed even from the figures. For / < 5 the agreement is still within an error of order 10"'^. 



More comments and examples can be found in |27, 28 
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4.6.3 Holes and close roots 

Let us now extend our investigation to situations with complex roots and consider the two particle 
states in more detail. Forgetting for the moment the breathers, we have to consider the two- 
soliton states. The soliton-antisoliton come in doublets so there are four different polarizations 
for two particle states of which the ss and ss have topological charge Q = ±2 instead ss and ss 
have zero topological charge. The first two states are expected to have exactly the same scaling 
function for energy and momentum because the sG/mTh theory is charge conjugation invariant. 
Instead there are two different situations for the neutral ss state, which have spatially symmetric 
and antisymmetric wavefunctions (denoted by (ss)+ and (ss)_, respectively). To separate the 
symmetric and antisymmetric part one simply has to diagonalize the 4x4 SG two particle S- 
matrix and see that it has 2 coinciding eigenvalues (equal to e^^^^\ corresponding to ss and ss), 
and two different eigenvalues in the Q = channel. 

Now we proceed to demonstrate that the IR limit restricts the possible quantum numbers of 
the complex roots. To simplify matters we consider only the repulsive regime p > 1. 

In the repulsive regime p > 1, the scaling function (ss)_ is realized as the solution to the NLIE 
with two holes (at positions /ii^2) and a complex pair at the position p±ia . In the IR limit we 
have 

Z(-i?) = I sinh-d + xi^- hi) +x{^ - h2) -xi"^- p-icr) -X{^- P + icr), 

Z{hl,2) = 2^/1,2 , ^^-^^^ 
Z{p ± ia) = 2itI^ . 

The quantization condition for the complex roots explicitly reads (we write down the equation 
only for the upper member of the complex pair, since the other one is similar) 



lsmh{p + ia) +x{p + icr - hi) +x{p + icr - h2) - x(2icr) = 27r/+ (4.36) 

Now observe that as / oo, the first term on the left hand side acquires a large imaginary part, 
but the right hand side is strictly real. The imaginary contribution should be cancelled by some 
other term. The function x is bounded everywhere except for isolated logarithmic singularities on 
the imaginary axis. For the cancellation of the imaginary part the argument 2ia of the last term 
on the left hand side should approach one of these singularities (similarly to the analysis in TEA 

ll)- 

In the repulsive regime, taking into account that for a close pair a < tt, the only possible 

vr 

choice for a is to approach ^ as / — > oo. The soliton-soliton scattering amplitude has a simple 

zero at = ivr with a derivative which we denote by C (the exact value does not matter). To 
leading order in I, the cancellation of the imaginary part reads 

icoshp + 3f?elog(7(^fT- I) =0, (4.37) 



from which we deduce 

exp(-/cosh(/))) , (4.38) 



TT 

^"2 



so the imaginary part of the complex pair approaches its infrared limit exponentially fast. This 
approach is modified by taking into account the finite imaginary contributions coming from the 
source terms of the holes and from the convolution term. These contributions lead to corrections 
of the order e~' and so they modify only the value of the constant C. 
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The behaviour of the real part near the singularity is: 

^Hexiz-in) ^0 if Kez < 0. 

Assume, for the moment, that a moves toward | from below (i.e. cr — ^ < 0). 
For the real part we get, again to the leading order 



^ex {p + i^- hi) + ^ex (p + - ^2) = 27r/, 
It can be shown that (in the repulsive regime p > 1) 

1 / TT 

1? 



c . 



(4.39) 



vr 



sinh — (i — 

2 log 



sinh — (i \- 'd 

p \ 2 



(4.40) 



where the the fundamental branch of the logarithm is taken into account. ^ is an odd monotonous 
function bounded by 



ie(^)i < ic(oo) 



7r(p - 1) 
2p 



(4.41) 



from below and above. This means that for any allowed value of 1^ the real position of the 
complex pair is determined uniquely and that 



l^^l< 



p-1 



2p 



(4.42) 



and since in the repulsive regime p > 1, the only possible choice is I(j = 0. Then the solution for 
p is 



hi + /i2 



(4.43) 



so it approaches the central position between the two holes (the corrections to this asymptotic 
are also exponentially small for large I). In fact, for the symmetric hole configuration Ii = — 12 
we expect /ii = — /i2 and p = to be valid even for finite values /. 



vr 



However, the above derivation is valid only if cr < — and so we do not cross the boundary of 

vr 

the analyticity strip of the x(^) function which is at 9m = vr. If cr > — , we must use a consider 
that: 



SRe x(-z — ^7r) 



vr 



if 



3f?ez > 0. 



The conclusion that a approaches ^ exponentially fast as Z — > 00 is unaffected by this change, 
but now the real part of x{2icr) is not zero but instead ±vr (modulo 27r). We choose the branches 
of the logarithm in such a way that Kex(i2icr) = =Fvr (the structure of the cuts is displayed in 
figure 4^). In this way we obtain that 



^2' 



and the asymptotic value of p is just as before. 
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Figure 4.6: The analytic structure of tlie counting function if cr > ^. The logarithmic cuts are indicated 
with the wiggly lines. The one lying in the upper half plane originates from the root in the lower half 
plane and vice versa. We also indicated the value of the discontinuity across the cuts. 



The S-matrix of the (ss)- configuration can be obtained by substituting the asymptotic values 

hi + h2 .vr , ,x 

p = ^ ' ^ = ^ 2 ^ ^ ^ 

into the expression for Z{'d) and considering now the quantization rules of the holes. We obtain 

Z{hi) = /sinh(/ii)-ilog5_(/ii-/i2) = 27r/i , (4.45) 

and a similar equation for the second hole, with 

5_(^) = JJ^^XX^^)^ (4-46) 

which is the correct scattering amplitude for the antisymmetric configuration of the soliton- 
antisoliton system. The amplitude has poles at i? = m{l — 2kp) , k = 1,2,... , corresponding 
to breathers with mass 

m2k = sin kpir . 



The equation ( [4.45 ) describes an approximation to the full NLIE valid for large I which is called 



the dressed or asymptotic Bethe Ansatz. There are three levels of approximations to the scaling 
functions: the full NLIE (which is in fact exact), the higher level Bethe Ansatz (HLBA) obtained 
from the NLIE by dropping the convolution term and the asymptotic BA. The difference between 
the BA and the HLBA is that while the former keeps the complex pair in its asymptotic position, 
the latter takes into account the corrections coming from the fact that a changes with / (to first 
order as given in (Ul) ). However, since the convolution term is of the same order in I as the 
dependence of a derived from the HLBA, the HLBA is not a self-consistent approximation scheme. 
Therefore we are left with the exact NLIE, which is valid for all scales and with the BA as its IR 
asymptotic form. 

The above conclusions about the (ss)_ state an be extended into the attractive regime as 
long as p > ^. At the point p = ^, however, the pair situated at the asymptotic position 
hits the boundary of the fundamental analyticity strip situated at inp. This phenomenon has a 
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physical origin: this point is the threshold for the second breather bound state which is the first 
pole entering the physical strip in the (ss)_ channel. To continue our state further, it requires to 
go to configurations having an array of complex roots of the first type (see ||8|]). Without going 
into further details, let us mention only that any such array contains a close pair plus some wide 
pairs depending on the value of p. 

The symmetric state (ss)+ can be obtained from a configuration with two holes and one 
selfconjugate wide root. Considerations similar to the above lead to the scattering amplitude 



cosh 



i) + in 
2p 



cosh 



•d — ITT 

2p 



which agrees with the prediction from the exact S-matrix. The poles are now at ■(? = Z7r(l 
l)p) , k = 0,1, . . . and correspond to breathers with mass 



(4.47) 



{2k + 



n^2k+i = sin 



{2k + l)p7r 



This configuration extends down to p = 1, where the selfconjugate root collides with the boundary 
of the fundamental analyticity strip at itt. This is exactly the threshold for the first breather which 
is the lowest bound state in the (ss)+ channel. For p < 1 the (ss)+ state contains a degenerate 
array of the first kind (see |^), which consists of a close pair, a selfconjugate root and some wide 
pairs, again depending on the value of p. 

We remark that it is easy to extend the above considerations to a state with four holes and 
a complex pair. In this case the only essential modification to the above conclusions is that the 
limit for the complex pair quantum number becomes 



p — 1 



when a < —, and 
2' 



< 



p 



p — 1 



(4.48) 



(4.49) 



TT 



when ^ > "^j since now there are four x sources coming from the holes. In the range 1 < p < 2 
this gives the same constraints 



0, 



a < 



2' 



a > 



TT 

2 

IT 



as for the (ss)_ state (for p > 2 it allows for more solutions, which however we will not need in 
the sequel). 

From the UV point of view, the results shown in this paragraph about the quantum numbers 
of the close roots are quite interesting. Consider the conformal dimensions ( [4.26| ) or the specific 
form ( 4.27| ). There is no proof that the terms N± give always the correct descendent numbers. In 



particular there is no proof that the sum of quantum numbers (/^ — 21^ — — is bounded 
from below (and nonnegative) . Then we don't know if the so called "vacuum" is really the hamil- 
tonian ground state. The IR analysis done until now shows exactly that states with only holes 
and with holes and close roots have the correct hamiltonian behaviour. 
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We now turn ourselves to the UV limits of the previous states. In particular we consider the 
state of two holes and a close pair (in the repulsive regime), which describes the antisymmetric 
soliton-antisoliton two-particle state. According to the results from the IR asymptotic, the quan- 
tum numbers of the complex roots can be or ±i. Therefore, we are left two possibilities: either 
the holes are quantised with integers or with half-integers. 

Let us start with the half-integer choice 6 = and suppose that one of the holes is right moving 
(which is the case if its quantum number Ii = /+ is positive) and the other one is left moving (i.e. 



I2 = I- < 0). Using the general formula ( 4.17 ) and = |, 5 = 0, by a simple calculation we 
obtain 

Q+(-oo) = Z+(-oo) = -27r^ , 

p+l 

which is valid for 1 < p < 4 (the complex pair remains central) . The other plateau values follow 
by oddity of the function Z . For the other half of the repulsive regime we have to include two 
new normal holes, one moving to the left and the other to the right, and two special holes that 
remain central (a justification for this will be given shortly) so that = 1. We then find 

Q+(-oo) = 47r^ , Z+(-oo) = -27r^ , 
p+l p+l 

which is valid for p > 4. In both cases the result for the UV weights turns out to be 
= +I±-- = +I±-- = Af±i 0) ± /± - - . 

The primary state is obtained for the minimal choice /+ = — /_ = | and coincides with a linear 
combination of the vertex operators g) which is correct for the state to be included in the 
spectrum of the sG/mTh theory and agrees with the behaviour of the {ss)_ state observed from 
TCS. If both of the holes move to the right or to the left, we obtain descendents of the identity 
operator, i.e. states in the vacuum module of the UV CFT. 



Figure |4/7| shows how the special holes are generated in this case. Starting from p < 4 and 
increasing the value of p the UV asymptotic form of the counting function varies analytically. 
Since the real roots/holes are quantised by half-integers, they are at positions where the function 
Z crosses a value of an odd multiple of it. As the plots demonstrate, the behaviour of Z is in 
fact continuous at the boundary p = 4: it is our interpretation in terms of the sources that 
changes, exactly because we try to keep the logarithm in the integral term of the NLIE in its 
fundamental branch. The price we pay is that we have to introduce two new normal holes (one 
moves left and the other one moves right) and two special holes which are central. Let us comment 
on the extension to the attractive regime. The complex root configuration changes only by the 
possible presence of wide pairs, which however have no effect on the plateau values ( 4.17D as their 
contribution can be absorbed in redefinition of the values k±. The other effect the wide roots 



have is to shift the terms S± by multiples of 2tt, but this affects only the integers N± in (|4.29| ). 
Therefore we again obtain states which are descendents of V(-|-i q) . 

For the state with integer quantization of the holes, we only give the result. If one of the holes 
is moving left and the other one to the right, we obtain 

i.e. some linear combinations of descendents of the vertex operators V(_|_i/2,o)- If both of the holes 
move to the left or to the right, we obtain other descendants of the same primary states. Again, 
this excludes the integer quantised states from the spectrum of sG/mTh theory. 
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4.6.3.1 Two holes and a selfconjugate complex root 

Taking the integer quantised state, the plateau values are given by the same formulas as for the 
case with the close pair above, since the plateau equation is identical. The only difference is that 
the numbers ( fl.lOl) take a nonzero value 

''W ~ ■ 

If the one of the holes is a right mover, while the other one is a left mover, the conformal weights 
turn out to be 

A± = ± /± - 1 = A(±i,o) ± /± - 1. 

If the two holes move in the same direction, we again obtain secondaries of the vacuum state. 

Concerning the extension to the attractive regime, one can again check that the plateau 
solution remains unchanged for the root configuration of two holes and two close roots, therefore 
the conformal family to which the state belongs remains the same, similarly to the case of the 
(ss)_ state. 

One can ask what happens if we quantize with half-integers. The result is, similarly to the 
case with the complex pair, that we obtain descendants of the operators V(±i/2,o)- Therefore 
we conclude that in this case the integer quantised configuration must be accepted, while the 
half-integer one is ruled out, in full accordance with the rule ( 4.31J ). 



4.6.4 Breather S'-matrices and IR limit 

In the infrared limit / — > oo the term I sinh(i?) develops a large imaginary part in the first deter- 
mination away from the real axis, forcing the close complex roots to fall into special configurations 
called arrays (we use the terminology of |^]). An array is a set of complex roots in which the 
roots are placed at specific intervals in the imaginary direction and have the same real part. In 
the attractive regime / sinh(T?)// is nonzero and so this is true for nonselfconjugate wide pairs as 
well (in the repulsive case wide roots do not have such driving force), while self-conjugate roots 
have a fixed imaginary part anyway. The deviation of the complex roots from their positions in 
the array decays exponentially with / (section |4.6.3| ). 

For the rest of this subsection, whenever it is not explicitly stated, we restrict ourselves to the 
attractive regime p < 1. The possible arrays fall into two classes: 



1. Arrays of the first kind are the ones containing close roots, which describe the polarization 
states of solitons. 
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There are two degenerate cases: odd degenerate arrays, which have a self- conjugate root at 

.7r(p + 1) 



^0 = + i- 



and accompanying complex pairs at 



7r(l - (2A; + l)p) , 
= ± i— ^- — , fc = 0, . . . 



1 

2p 



and even degenerate ones, which only contain complex pairs, at the positions 



1 ■ 
2p. 



These arrays always contain exactly one close pair. The odd degenerate arrays in the 
repulsive regime reduce to single self- conjugate roots and the even degenerate ones to a 
single close complex pair. 

2. Arrays of the second kind describe breather degrees of freedom. The odd ones contain a 
self- conjugate root 



and wide pairs as follows: 
where 



7r(l - (2A: + l)p) 



k = 0, . . . , s , 



0<s< 



while the even ones only contain wide pairs 



1 

2P. 



9mi?fc = i?±z^^LM ,k = 0,...,s, 

and s runs in the same range. They correspond to the (2s -|- l)-th breather ^25+1 and the 
(2s -|- 2)-th breather i?2s+2, respectively. 

As one can see, arrays of the second kind become degenerate ones of the first kind, if we analytically 
continue increasing p. The reason is that breathers are of course soliton-antisoliton bound states, 
while degenerate arrays of the first kind describe scattering states of a soliton and antisoliton, as 
we will see shortly. 

One can compute the energy and momentum contribution of a array of the second kind cor- 
responding to the breather Bs- The energy-momentum contribution turns out to be 



TTST) 

2M sin (cosh 1?, sinh ??) , 



(4.50) 



where is the common real part of the roots composing the array. This is just the contribution 
of a breather Bg moving with rapidity i?. Arrays of the first kind do not contribute to the energy- 
momentum in the infrared limit, which lends support to their interpretation as polarization states 
of solitons. 
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Now we proceed to show that with the above interpretation the NLIE correctly reproduces 
the two-body scattering matrices of sine-Gordon theory including breathers. 

Let us start with breather-soliton matrices. The Bethe quantization conditions for a state 
containing a soliton (i.e. a hole) with rapidity t?i and a breather Bs with rapidity i?2 take the 
following form in the infrared limit. For the hole we get 

s 

Z{{)i) = M sinhi?i - ^ xiii^i -^2- ipk) = 27r/i, 

A:=0 

where we denoted the prescribed imaginary parts of the roots in the array B^. by pk. Here we used 
x(0) = to eliminate the source term for the hole. Now we can compute the second determination 
of X as in ( PI) . Now it is a matter of elementary algebra to arrive at 



Z{^i) = M sinh {>i - i log SsBs {^1 - ^2) = 27r/i , 

where SsBsi'^^i ~ ^2) is the soliton-breather 5-matrix conjectured in 

One can start with the breather quantization conditions, too. Writing 

Z {^2 + ipk) = M sinh {d2 + ipk)n + Xi02 -^1+ ipk) + ■■■ = 2vr/f ^ , A; = 0, . . . , s (4.51) 

(the dots are terms due to wide root sources themselves, which cancel out up to multiples of 27r 
in the next step). Summing up these equations one arrives at 

27W sin ^sinh(T?2)-« log 55B,(??2-^9i) = 2^/2 , 
where I2 is essentially minus the sum of the quantum numbers of the wide roots composing the 



array (shifted by some integer coming from summing up the terms omitted in equation ( 4.51 )). 

Using a similar line of argument we also reproduced the breather-breather 5-matrices by 
writing down the Bethe quantization conditions for a state with two degenerate strings Bg and 
Br of the second kind. One has to be careful that when Z{'d) contains wide root sources which 
are expressed in terms of Xiii^)} the second determinations of these terms will appear in Z('i?)for 
large Qm-d, i.e. terms that can be written roughly like {xil{'&))ii, as in ( 3.39| ). 



Scattering state of a soliton and an antisoliton can be described by taking two holes and a 
degenerate array of the first kind. There are two possibilities now, corresponding to scattering in 
the parity-odd and parity-even channels. Following the procedure outlined in section 4.6.3, we were 



once again able to reproduce the corresponding scattering amplitudes. The results presented here 



together with those of section 4.6.3 exhaust all two-particle scattering amplitudes of sine-Gordon 



theory, both in the repulsive and attractive regime. 



4.6.5 Some examples of breather states 

So we proceed to take a look at the simplest neutral excited state, which is the one containing a 
first breather Bi at rest. The source term to be written into the NLIE turns out to be 

cos I smh[v — wr) 

g{i^ - wr) = -xiii-^ - w) = -Hog , (4.52) 

cos ^ + ^ sinh(T? — wr) 

where wr is the real part of the position of the self-conjugate root (the imaginary part is wj = 
tt{p+ l)/2. In this case wr = since the breather has zero momentum. There is no need to look 
at the Bethe quantization condition as the root does not move due to the left-right symmetry 
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I 


TCS 


NLIE 


TCS 


NLIE 


1.0 


2.38459 


n/a 


1.84996 


n/a 


1.5 


1.68168 


n/a 


1.30438 


n/a 


2.0 


1.35391 


n/a 


1.05038 


n/a 


2.5 


1.17420 


n/a 


0.91152 


n/a 


3.0 


1.06692 


1.0655810539 


0.82903 


0.8285879853 


3.5 


0.99985 


0.9980153379 


0.77773 


0.7771857432 


4.0 


0.95664 


0.9542867454 


0.74499 


0.7443106400 


4.5 


0.92845 


0.9254766107 


0.72381 


0.7229895177 


5.0 


0.91000 


0.9063029141 


0.71006 


0.7090838262 



Table 4.2: The first breather state at p = - and P — g- Energies and distances are measured in units of 
the soliton mass M, and we have subtracted the predicted bulk energy term from the TCS data. 



of the problem. This state is quantized with integer Bethe quantum numbers, i.e. S 



1, as in 



Table ^ presents the energy values obtained by iterating the NLIE in comparison to results 
coming from TCS, at the values p = 2/1 and p = 2/9. 

The table shows that iteration of the NLIE fails for values of I less than 3 (the actual limiting 
value is around 2.5). What is the reason? 

We plot the counting function Z on the real line for / = 3 and Z = 2 in figure 4^. In a first 
approximation we can safely neglect the integral term for these values of the volume to see the 
qualitative features that we are interested in. What we see is that the behaviour of the function 




.15 




0.1 




.05 




-1 y6.5 


\^ n/s 1 


/ -0.05 




/ -0.1 




/ -0.15 




Z{d) at Z = 3 and Z = 


= 2, respectively. 



changes: its derivative changes sign at the origin. As a result, two new holes appear where the 

new real zeros of the function are. But the topological charge remains zero, due to the fact that 

now we have a special root at the origin and so Ns = 1 and Nh = 2. The two new holes do not 

give us any new dynamical degrees of freedom: their quantum numbers are fixed to be and so 

their positions are uniquely determined. 

Of course when we calculated the UV dimension, the appearance of the new sources had to 

1 

be taken into account. It turns out that for - < p < 1 the two holes are left/right movers, while 
for p < — they remain central. Calculating the UV conformal dimensions we obtain 



A 



± 



P 

P + l 



(4.53) 



so the ultraviolet limit of this state is a linear combination of the vertex operators V±i_o of the 
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c = 1 UV CFT. This is in perfect agreement with the TCS calculations performed by us as in 



table 4.2 



How does this change of sign in the derivative of Z affect the iteration scheme for the NLIE 



(see section [4.4|) ? The two new zeros of Z{'d), which is a complex analytic function apart from 
logarithmic branch cuts, actually correspond to singularities of the logarithmic term in the NLIE 
( ^.54 ). They come along the imaginary axis in the plane as we decrease /, and at a certain 



point they cross our integration contour which runs parallel to the real axis at distance tj. As 



they make the logarithmic term in our NLIE (|3.54| ) singular, they blow up our iteration scheme. 



After reaching the origin of the 1) plane (at exactly the radius where the derivative of Z becomes 
0), they continue to move along the real axis, which corresponds to crossing a square root branch 
cut. 

This problem of numerically solve NLIE appears in all the cases where special roots/holes are 
taken into account (e.g. the "-" vacuum, the two holes-close quantized with "-" and so on). 

We do not go into details here as this problem is currently under investigation^. We just 
remark that these issues prove to be highly nontrivial and for the time being, unfortunately, they 
prevent us from having a reliable numeric scheme for the NLIE below the critical volume. 

One can estimate the volume where the slope of the counting function changes sign by ne- 
glecting the integral of the logarithmic term. The result is 

2 

^critical p^f (4.54) 

COS- 

which gives a value of around 2.22 for p = 2/7 and 2.13 for p = 2/9. The actual limiting value is 
a bit higher, partly due to the finite value of t] used in the iteration program and partly because 
the iteration already destabilizes when the singularities come close enough to the contours. It 
must also be noted that the integral term cannot eventually be neglected when the singularities 
are close to the contour, which is an additional reason why ( 4.54| ) is just a crude estimate. 



We make a short digression to examine the UV limit of the second breather i?2- The second 
breather at rest is described by a wide pair at positions 

-d = ±i- . 
2 

Calculating the UV conformal dimension we get 

p 



A+ = A" 



P + l 



which turns out to be the same as that of the first one ( [1.53 ), i.e. this state must originate from 



the other linearly independent combination of the vertex operators V±ifl in the ultraviolet. This 



is again in perfect agreement with TCS and confirms a result by Pallua and Prester |g9[ who 
used XXZ chain in transverse magnetic field to regularize sine-Gordon theory. They calculated 
scaling functions numerically on a finite lattice for several concrete values of p, and arrived at this 
conclusion by looking at the numerical data. However, our method to compute UV dimensions 
gives us an exact analytic formula and therefore much stronger evidence. This result is interesting 



because it invalidates a conjecture made previously by Klassen and Melzer |20| who identified the 
second breather linear combination of V±2,o- 

To close this section, we present the lowest lying example of a two-breather state, containing 
two Bi particles with zero total momentum. It turns out that this is a state for which the numerical 
iteration of the NLIE is not plagued with the problem found above for the first breather. Locality 

^ Work in progress in collaboration with P. E. Dorey and C. Dunning, Durham. 
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constrains the state to be quantized with half-integers and for lowest energy the quantum numbers 
of the self-conjugate roots must take the values 



We remark that in contrast to the case of holes, the self-conjugate root with I > moves to 
the left, while the one with / < moves to the right. This is due to the fact that the second 
determination of Z is in general a monotonically decreasing function on the self- conjugate line. In 
order to determine the position of the two self- conjugate roots we need the second determination 
of Z. The second determination of the self- conjugate root source turns out to be 



i sin np — sinh i) 
i sin np + sinh i) 



Up to some signs, this is just the phase shift which arises when two breathers scatter on each 
other, which is exactly why the IR analysis gives the correct scattering amplitude. Using this 



formula, we obtained the numerical data presented in table 4.3. The UV limit for this state can 



/ 


TCS 


NLIE 


1.0 


12.1601 


12.159257 


1.5 


8.20139 


8.2006130 


2.0 


6.24771 


6.2465898 


2.5 


5.09489 


5.0937037 


3.0 


4.34132 


4.3397021 


3.5 


3.81513 


3.8129798 


4.0 


3.43020 


3.4275967 


4.5 


3.13912 


3.1357441 


5.0 


2.91308 


2.9089439 



Table 4.3: The two-breather state at p 



be calculated from NLIE to be a symmetric first level descendent of the vacuum with weights 

A+ = A- = 1 , 

which agrees with TCS. (Note that this descendent exists due to the fact that there is a U{1)l x 
U{1)r Kac-Moody symmetry at c = 1: this state exactly corresponds to the combination of the 
left and right moving currents J J.) 



4.7 a twist and minimal models (the ground state) 



It is known from ||23|| that by twisting the Bethe Ansatz Equations of the six vertex model as 



indicated in section |2^, the twisted model shows the critical behaviour of conformal minimal 
models (the untwisted critical behaviour corresponds to c = 1). Moreover, in recent papers, Al. 



Zamolodchikov has put forward the idea of modifying sine-Gordon theory by a twist a |30| to deal 
with conformal minimal models. As a consequence of this two ideas, the Bethe equations used 
to obtain NLIE contained the twist ( ^.37 ) to. In the NLIE, as shown in section |3.5| , the twist is 
parametrized by a ( 3.42| ). 

Looking at the ground state, that in analogy with sine-Gordon is expected to be a sea of real 
roots, the source in NLIE is put to zero §{'&) = ( |337| ). 
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In analogy with the sine-Gordon ground state, we can choose half-integer quantization rule 
with 6 = 0. The expression for conformal dimensions ( [4.26|) gives: 

24 Ap+1 Vvr/ 



corresponding to an effective central charge 



(the effective central charge is defined in ( |4.5| )). Only in the unitary models the c is the Virasoro 
central charge. Furthermore, it is well-known that the perturbation of the Virasoro minimal 
model Vir{r,s) by its relevant primary operator <&(i,3) is integrable and is described by an RSOS 



restriction of sine-Gordon theory |31| with 

p=—!—. (4.56) 

s — r 

We will use for this model the shorthand notation Vir{r,s) + '^(1,3). Using the rule suggested 
in (|3.60| ) for u) and the expression ( ^.42 ) for a gives: a = vr/r (the additional integer terms are 
excluded, for the moment). Then 

6 

c=l-- , 
rs 

which is exactly the effective central charge of the minimal model Vir{r,s). Therefore one can 
expect that the twisted equation describes the ground state of the model Vir{r,s) + ^(1,3)- In 
fact, Fioravanti et al. ||3^ calculated these scaling functions for the unitary case s = r + 1 and 
showed that they match perfectly with the TEA predictions already available. Moreover, choosing 
the following values for the twist 

a = ±— , fc = l...r-l (4.57) 

r 

they obtained the conformal weights of the operators ^{k,k) > k = 1 . . .r — 1 in the UV limit (the 
sign choice is just a matter of convention). In our notation, ^{q,q') denotes the primary field with 
conformal weights 

h^ = h-= . (4.58) 

4sr ^ ' 

The models Vir{r,s) + $(1,3) have exactly r — 1 ground states. In fact, one can see from the 
fusion rules that the matrix of the operator $(1,3) is block diagonal with exactly r — 1 blocks in 
the Hilbert space made up of states with the same left and right primary weights. In each of these 
blocks, there is exactly one ground state and for the unitary series s = r-|- 1, it was conjectured in 
1^] that their UV limits are the states corresponding to ^(^k,k)- One can check that in the general 
nonunitary case the twists ( 4.57| ) correspond in the UV limit to the lowest dimension operators 



among each of the r — 1 different blocks of primaries (see explicit examples later). 

These ground states are degenerate in infinite volume, but for finite I they split; their gaps 
decay exponentially as / — > 00. In the unitary case, they were first analyzed in the context of the 
NLIE in where it was shown that the NLIE predictions perfectly match with the TEA results 
already available for the unitary series. 

However, ground states for nonunitary models have not been treated so far and therefore now 
we proceed to give examples of that. The models we select are the ones that will be used for 
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comparison in the case of excited states as well. The first is for the scaling Lee- Yang model 



Vir{2, 5) + ^(1,3), for which we have also given data from TCS |15| and TBA |16| for comparison 
(table 0). 

Mb = 2M sin ^ = V2,M 

is the mass of the fundamental particle of the Lee- Yang model (this is more natural here than 
using the mass M. of the soliton of the unrestricted sine-Gordon model as a scale, since the soliton 
disappears entirely from the spectrum after RSOS restriction) and call 

Ib = MbL = Vsi 

where I is the variable appearing in NLIE. There is only one independent value of the twist, which 
we choose to be 

vr 

Here and in all other subsequent calculations the TCS data were normalized using the analogue 



of the coupling-mass gap relation ( 1.18 ) from |25]. There is only one ground state in this model. 



Ib 


TCS 


NLIE 


TBA 


0.1 


-2.0835015786 


-2.0835015787 


-2.0835015786 


0.5 


-0.3803475256 


-0.3803475281 


-0.3803475281 


1.0 


-0.1532068463 


-0.1532068801 


-0.1532068801 


1.5 


-0.0763483319 


-0.0763484842 


-0.0763484842 


2.0 


-0.0406269362 


-0.0406273676 


-0.0406273676 


2.5 


-0.0222292932 


-0.0222302407 


-0.0222302407 


3.0 


-0.0123492438 


-0.0123510173 


-0.0123510173 


3.5 


-0.0069309029 


-0.0069338817 


-0.0069338817 


4.0 


-0.0039198117 


-0.0039244430 


-0.0039244430 


5.0 


-0.0012721417 


-0.0012816882 


-0.0012816882 



Table 4.4: The vacuum of the Virasoro minimal model Vir{2,5) perturbed by $(1,3). The energy and 
the volume are normalized to the mass of the lowest excitation Mb, which is the first breather of the 
unrestricted sine-Gordon model. The TCS data shown have the predicted bulk energy term subtracted. 



which corresponds to the primary field with conformal weights 

A+ = A" 



1 

'5 ' 



which is in agreement with TBA and TCS predictions. We have also found a perfect agreement 
for the models Vir{2,7) + ^(1,3) and Vir{2,9) + '^(1,3), but we do not present those data here. 
We remark that the TCS for the minimal models converges much better than the one for c = 1 
theories: all TCS data in table |4j and subsequent ones were produced by taking a few hundred 
states and in some fortunate cases (e.g. the ground state of the scaling Lee- Yang model for 
small values of I) we were able to produce data with up to 9 — 10 digits of accuracy! The better 
convergence meant that all the computation could be done with the computer algebra program 
Mathematica, greatly simplifying the programming work. 

All the models of the class Vir{2, 2n + 1) -|-'&(i^3) have only one ground state. For models with 
two ground states, we can take a look at yir(3,5) (Vir{2>,7) was taken into account in p^]). For 
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Vir{3,5) the ultraviolet spectrum is defined by the following Kac table, where the weight ( 4.58 ) 
of the field ^(^k,l) is found in the k-th row and l-th column. 






i 

20 


i 

5 


3 
4 


3 
4 


i 

5 


i 

20 






The two blocks of the perturbing operator <&(i,3) are defined by the fields {$(1,2) > ^(1,4)} 
{$(11) , $(13)}, respectively. The ground states correspond in the UV to the operators $(1,2) and 
^(1,1), as can be checked directly using formulae ( W^) , (HI) and 

We also have TBA data to compare with, using the TEA equation written by Christe and 



Martins |35]. The lower-lying ground state is obtained directly from their TBA, while for the 
other we used the idea of Fendley of twisting the TBA equation ||36|. The numerical results are 



presented in tables and 4.6 



In the case Vir{3,7) treated in |34|, it was possible only to have a comparison with TCS 
results, but it still looked pretty convincing. To summarize, we now have sufficient evidence 



I 


TCS 


NLIE 


TBA 


0.1 


-3.074916 


-3.0749130189 


-3.0749130190 


0.3 


-0.944161 


-0.9441276204 


-0.9441276204 


0.5 


-0.509764 


-0.5096602194 


-0.5096602194 


0.8 


-0.265436 


-0.2651431026 


-0.2651431026 


1.0 


-0.186038 


-0.1855606546 


-0.1855606546 


1.5 


-0.087300 


-0.0861426792 


-0.0861426792 


2.0 


-0.045910 


-0.0437473815 


-0.0437473815 


2.5 


-0.026746 


-0.0232421927 


-0.0232421927 


3.0 


-0.017868 


-0.0126823057 


-0.0126823057 


4.0 


-0.013546 


-0.0039607326 


-0.0039607326 



Table 4.5: One of the two ground states of the Virasoro minimal model Vir{3,5) perturbed by $(1,3), 
corresponding to a = ^. The energy and the volume are normalized to the mass of the lowest excitation, 
which is the soliton of the unrestricted sine-Gordon model. The TCS data shown have the predicted bulk 
energy term subtracted. 



I 


TCS 


NLIE 


TBA 


0.1 


3.117844 


3.1178476855 


3.1178476853 


0.3 


0.985360 


0.9853990810 


0.9853990810 


0.5 


0.540427 


0.5405470784 


0.5405470784 


0.8 


0.282725 


0.2830552991 


0.2830552991 


1.0 


0.197143 


0.1976769278 


0.1976769278 


1.5 


0.089277 


0.0905539780 


0.0905539780 


2.0 


0.042960 


0.0453290013 


0.0453290013 


2.5 


0.019978 


0.0238075022 


0.0238075022 


3.0 


0.007209 


0.0128843786 


0.0128843786 


4.0 


-0.006592 


0.0039866371 


0.0039866371 



Table 4.6: The other ground state of the Virasoro minimal model Vir{3,5) perturbed by ^(1,3), cor- 
responding to a = The energy and the volume are normalized to the mass of the lowest excitation, 
which is the soliton of the unrestricted sine-Gordon model. The TCS data shown have the predicted bulk 
energy term subtracted. 
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to believe that the a-twisted NLIE describes the correct scaling functions for ground states of 
minimal models perturbed by ^(1,3) in unitary and nonunitary case. However, the NLIE for sine- 
Gordon is known to work for excited states as well. But how do we get the excited state spectrum 
of the minimal models now? 



4.8 a twist and minimal models: excited states 
4.8.1 The choice of a 

From now on we restrict ourselves to the case of neutral (i.e. S = 0) states. It is easy to see that 
even for states with a zero charge the relation between a and to is highly nontrivial. 
Now we proceed to show that choosing the value of u! as 



p+l 



where k is integer, we can reproduce all the required values of a listed in equation ( [4.57|) . Observe 
that this expression for ui is different from ( |3.60|), but totally equivalent, as will be clear soon. 
First of all, we substitute the value of p from (|4.5^) to obtain 



k{s — r)TT 



Since r and s are relative primes, the independent values of uj mod vr can be written as 



CO 



Itt , 

— , / = 0,... . 

s 



For = 0, we can rewrite the formula (|3.42[) as follows 



a 



Itt 2r — s 

\ TT 

r 2r 



1 r 

2+i 



1 / 

2 ~ s 



We are interested only in the value of a mod vr, since using the parameter 6 one can effectively 
shift a by TT. This leaves us with the formula 



a 



In sir 
r 2r 



1 r 

2 + 1 



1 / 

2 ~ s 



The first possibility is that ^ < f , which simply gives us the values 



a 



Itt 
r 



When s is even, we can have ^ = §, which gives us a = 0. Finally, when / > |, we get the values 



a 



(l - s)7r 



It is easy to check that these formulae reproduce every value 



niT 

a = — mod vr 

r 



at least once, as required by ( 4.57| ), using the fact that s > r and that the values above form an 
uninterrupted sequence of s numbers (or when s is even, of s — 1 numbers, the zero repeated) with 
equal distances ^. 
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As we have already seen in the previous section, all the values 

fcvr , , , 

a = — , k = l,... ,r-l (4.59) 
r 

are necessary to reproduce correctly the r — 1 ground states of the model Vir{r, s) + ^*(i,3)- 

The twisted lattice Bethe Ansatz was analyzed by de Vega and Giacomini in ||3^. On the 
lattice, passing from the sine-Gordon model to the perturbed Virasoro model amounts to going 



from the six-vertex model to a lattice RSOS model. In ||37|] it was shown that to obtain all the 



states of the RSOS model it is necessary to take all the twists 

kir 

ui = mod vr 

p+l 

into account. The fact that not all these twists correspond to inequivalent values of a and so to 
different physical states is a consequence of the RSOS truncation. 

To close this section we remark that the parameter a drops out of the second determination of 
Z in the attractive regime. This is important because as a consequence the IR asymptotics of the 
breather states does not depend on a and so the S-matrices involving breathers are unchanged. 
In fact, scattering amplitudes between solitons and breathers remain unchanged too, as can be 
seen from examining the argument that we used to derive them in section ^1.6.4 . This matches 



with the fact that the RSOS restriction from sine-Gordon theory to perturbed minimal models 



does not modify scattering amplitudes that involve two breathers or a breather and soliton |31] 



4.8.2 The UV limit 

There is in fact a very simple intuitive argument to show that the states we get from NLIE with 
an a twist are related to the minimal models in the UV limit. Let us first recall that the UV limit 
of the a = NLIE yields the vertex operators V(^n,m) their descendants ( [1.27 , 4.28 ). 



Let us look at the conformal dimensions. First note that because the value of a for a minimal 
model is never a multiple of vr, one does not expect central sources in the UV limit but only 
right/movers or left movers. Indeed, in all the examples of sine-Gordon with central sources, 
the left-right symmetry (i.e. Z{—'&) = —Ziif) on the real axis) of the NLIE was crucial. This 
symmetry, however, only holds for a = or vr. As a consequence we have 

5° = 0, 5 = 5+ + 5" , 

and in addition Q+(— oo) = Q_(-|-cx)), which means that only one-plateau systems are allowed in 
the twisted case. 



Using the formulas (|4.2q , [l.27| ) for the UV limit of the NLIE, one can see that introducing 
a 7^ is equivalent to shifting the quantum number n± to n± + . One has to be careful that 
since the value of the central charge is shifted from 1 to the one of the minimal model, we have to 
take this shift into account when computing the conformal weight from the leading UV behaviour 
of the energy level ( |4.5| ) . We put in ( [1.26D the value 

a = — , j = l,...,r-l, 

r 

as in ( [1.59| ), and the central charge 

6(r - sf 



c= 1 



rs 
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of the minimal model Vir{r,s), as in section 4.7. By a computation similar to the one done 
in section 4^ to obtain (|4.27| , [4.28 , 4.30 ) one can obtain the following expression for conformal 
dimensions: 



{Is - I'j^rf - (s - rf 

4rs 
/ = S 



(4.60) 



± 5 



+ 2k± 



2k^ 



2(5-2S±). 



The expression of N± is exactly the same that appears in ( [4.30 ). Moreover, from the formulae 
( [4.26 , 4.60| ) it is also clear that in general 



2fA^ 



2S- mod 1 
vr 



because only one plateau systems can take place, and so we see that general charged states will have 
fractional Lorentz spin. Actually, it is known that charged states in the models Vir{r,s) + <&(i,3) 
generally have fractional Lorentz spin p9| . 

Comparing the first line of ( 4.60 ) with ([4.58D one observes that it can represent the conformal 
dimensions of an operator ^{q,q') only if some conditions are verified. The most general one is 

(/s — ZV)^ = (gs — qr^ + 2rs ■ integer, l<q<r — 1, l<q'<s — 1 

(the /, I' are not integer numbers in general). Observe that it is an arithmetic (Diofantine) equation 
in q, q' and the unknown integer. 

It is convenient to treat some specific cases, because, as for the sine-Gordon case, there is no 
proof about the general behaviour of l'^ and A^±. 

In the unitary p = r, s = r + 1 and neutral case 3 = the resulting conformal weights take 
the form 



{2np + 



1 



- - 4p(p + i) ■ <«" 

We see that this is the weight of the field $(i,/_2n) in the minimal model Vir{p,p + 1), however, 
in order not to overflow the Kac table, the range of n must be restricted as 

l<l-2n<p . 

For charged states, a similar calculation can be performed. 

By an inspection of the formula (|4.26| ), for any state with 5 = 0, the property to have only 
one plateau implies that 

A+ - A~ 

is integer or half-integer. In fact, choosing the quantization rule and the parameter 6 in an 
appropriate way, one can ensure that this difference is integer (see IH, ^]). This means that 
the UV limit of any neutral state is either a field occurring in the ADE classification of modular 
invariants [38| or (in case we choose 5 so that — A~ is half-integer) it is a field from a fermionic 
version of the minimal model |20|. 

Until now it is not clear whether the weights actually stay inside the Kac table, for which in 
the unitary case one must require 

2k + 45+ 



1 < / 



< p . 



Due to the fact that the configuration of sources in the UV may be very non-trivially related to 
the one in the IR, this condition is very hard to check in general, but no concrete examples that 
we calculated have ever violated this bound. 



4.9. CONCRETE EXAMPLES OF EXCITED STATES 
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4.9 Concrete examples of excited states 

4.9.1 The Vir{2, 2n + 1) + $(1, 3) series 

Let us start with examining the scahng Lee- Yang model Fir (2, 5) + ^*(i,3)- There is only one 
independent value of the twist which we choose as 

vr 

2 ' 

since we have a single ground state, and as a result there are no kinks in the spectrum. We fix 
the value of a as above, so we still have a freedom of choosing 6. This can be done by matching 
to the UV dimensions: if for a certain state we choose the wrong value of 5, we find a conformal 
dimension that is not present in the Kac table of the model. 

The excited states are multi-particle states of the first breather of the corresponding unre- 
stricted sine-Gordon model, which has 

2 

P=3- 

Now one can calculate the state containing one particle at rest. We find the numerical data 



presented in table [1.7 . 

It turns out that as we decrease /, the self-conjugate root starts moving to the right. It does 
not remain in the middle like in the a = case, which is to be expected since for nonzero a we 
have no left /right symmetry. However, the total momentum of the state still remains zero due to 



a contribution from the integral term in momentum equation (|3.59D . One can see that once again 



h 


TCS 


NLIE 


0.1 


23.05277 


n/a 


0.5 


4.679779 


n/a 


1.0 


2.447376 


n/a 


1.5 


1.748874 


n/a 


2.0 


1.430883 


n/a 


2.6 


1.238051 


1.238012(#) 


3.0 


1.164321 


1.164319 


3.5 


1.105220 


1.105196 


4.0 


1.068256 


1.068237 


5.0 


1.029356 


1.029348 



Table 4.7: The first excited state of the scaling Lee- Yang model. The energy and the volume are 
normalized to the mass of the lowest excitation, which is the first breather of the unrestricted sine-Gordon 
model. The TCS data shown have the predicted bulk energy term subtracted. 

we have the phenomenon noticed in the case of the first breather of sine-Gordon theory, namely 
the appearance of the special root and its two accompanying holes, so the iteration breaks down 
again around I = 2.5. The (^) in the table 4.7 written after the NLIE result for / = 2.6 means 



that due to the fact that the singularities corresponding to the new holes and the special root 
are just about to cross the contour and upset the iteration scheme, the NLIE result becomes less 
precise. We will use this notation on later occasions too. In any case, the agreement still looks 
quite convincing. 

Let us now look at the UV spectrum of the model. We know that the Lee- Yang model contains 
only two primary fields, the identity I and the field ip with left/right conformal weights 

1 

~5 ■ 



A+ = A" 
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In fact, the ground state of the massive model corresponds to (/? in the UV limit. One can compute 
the UV limit of the first particle from the NLIE too, taking into account the appearance of the 
special root and the holes. It turns out that the special root and one of the holes moves to the 
left together with the self-conjugate root, while the other hole moves to the right. The result is 



A" 



, 



i.e. the identity operator I, which fits nicely with the TCS data (see also ||15| ). 

Let us look now at moving breathers. If the self-conjugate root has Bethe quantum number 
/ = 1, the corresponding state will have momentum quantum number 1, i.e. 

R ' 



and in the UV A+ — A^ = 1. One can note from the numerical data presented in table |4^ that 
the special root does not appear here. The reason is that the self- conjugate root moves to the left 
and the real part of its position i9 is given to leading order by 



sinh(3f?e??) 



2ttI 



B 



As a result, the contribution to the derivative of Z from the /sinhi? term remains finite when 
I ^ 0. In the previous example of the particle at rest the left-moving nature of the self-conjugate 
root when / = does not prevent the occurrence of the breakdown in the iteration scheme: since 
its Bethe quantum number is zero, it does not move fast enough to the left in order to balance 
the negative contribution to derivative of Z coming from the self- conjugate root source. At the 
moment we have no way of predicting analytically whether or not there will be specials in the UV 
limit: we just use the numerical results to establish the configuration for the evaluation of UV 
weights, supplemented with a study of the self-consistency of the solution of the plateau equation 
( 4.16|) . The UV dimensions for the moving breather turn out to correspond to the state L^iip. 



Ib 


TCS 


NLIE 


0.1 


60.75048 


60.74682 


0.5 


12.20618 


12.20561 


1.0 


6.182516 


6.182363 


1.5 


4.202938 


4.202915 


2.0 


3.231734 


3.231640 


2.5 


2.662186 


2.662110 


3.0 


2.292273 


2.292231 


3.5 


2.035552 


2.035530 


4.0 


1.848892 


1.848849 


5.0 


1.599792 


1.599762 



Table 4.8: The one-particle states with Lorentz spin 1 of the scaling Lee- Yang model. The energy and the 
volume are normalized to the mass of the lowest excitation, which is the first breather of the unrestricted 
sine-Gordon model. The TCS data shown have the predicted bulk energy term subtracted. 

One can similarly compute the UV dimensions for some other excited states. For example, 
the two-particle states with half-integer Bethe quantum numbers h > 0, I2 < for the two 
self-conjugate roots are found to have 



A+ 



A" 



h + 



1 
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in agreement with TCS data which show that they correspond in the UV to descendent states of 
if. The first such state with quantum numbers 



h 



corresponds in the UV to L^iL^iip and is given numerically in table 



/ 


TCS 


NLIE 


0.1 


123.583 


123.5693 


0.5 


24.7806 


24.77936 


1.0 


12.4870 


12.48635 


1.5 


8.42926 


8.428693 


2.0 


6.42931 


6.429201 


3.0 


4.48444 


4.484209 


4.0 


3.56367 


3.563519 


5.0 


3.04899 


3.048881 



Table 4.9: The lowest lying zero-momentum two-particle state in the scaling Lee- Yang model as computed 
from the NLIE and compared with TCS. 

The lowest lying three-particle state of zero momentum, with Bethe quantum numbers (—1, 0, 1) 
corresponds to the left/right symmetric second descendent of the identity field, i.e. to the field 
TT, where T denotes the energy-momentum tensor. This is very interesting, since from experience 
with NLIE UV calculations one would naively expect this to be a first descendent (descendent 
numbers are usually linked to the sum of Bethe quantum numbers of left/right moving particles 
and this state is the lowest possible descendent of the identity I). However, the field L_iZ_i][ is 
well-known to be a null field in any conformal field theory. 

The above correspondences are again confirmed by comparing to TCS (see the wonderful 
figures in fl^). In general, one can establish the rule that states with odd number of particles 
must be quantized by integers (J = 1), while those containing even number of particles must be 
quantized by half-integers {6 = 0) in order to reproduce correctly the spectrum of the scaling 
Lee- Yang model. 

We conducted similar studies for the models Vir{2, 7) + ^{1,3) and Vir{2, 9) -|- ^(1,3) and 
found similarly good agreement with TCS data. For the first one-particle state of the model 
Vir{2, 7) -|- *I'(i,3) we also checked our results against the TBA data in the numerical tables of |40| 
and found agreement with the TBA results. 

Given the choice of a above, the correct rule of quantization in all of the models Vir{2, 2n + 

1) + $(1,3) is 

5 = Msc mod 2 , 

where Msc is the number of self- conjugate roots in the source corresponding to the state. This is 
exactly the same rule as the one established for pure sine-Gordon theory in |^. In the presence 
of the twist, such a rule of course has meaning only together with a definite convention for the 
choice of a. 



4.9.2 One-breather states in the Vir{?i, 7) case 

It is interesting to note that in the case of Vir{?>,n) models, all the neutral states must come in 
two copies, since they can be built on top of either of the two ground states. We take the example 
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of the Vir{3,7) model and the states corresponding to a breather at rest. We have 

3 

but now there are two inequivalent values for the twist 

TT 27r 
" ~ 3 ' T ' 



When a = 0, we can calculate the critical value of I to be IcriUcai = 5.23 using (4.54). In this 
case, the twist helps a bit, because it makes the self-conjugate root a left mover; it is intuitively 
clear that the bigger the twist, the more it lowers the eventual value of IcrticaU which is in accord 



with the numerical results of table |4.10| . From the TCS data one can identify that breather #1 is 
really the one-particle state in the sector of ground state ^1 (a = vr/3), while breather #2 is in 
the sector over ground state #2 (a = 27r/3). 

A direct calculation of the conformal weights gives the following resultsiA"*" = A~ = ^for 
breather #1 and A"*" = A~ = Ofor breather #2, which are in complete agreement with the TCS 
data. We also checked the two different states containing two breathers with Bethe quantum 
numbersli = \ , I2 = —\ ,and found an equally excellent numerical agreement with TCS. Just 
like in the case of sine-Gordon and scaling Lee- Yang model, for these states one can continue the 
iteration of the NLIE down to any small value of I, although at the expense of a growing number 
of necessary iterations to achieve the prescribed precision. 





breather 1 


breather 2 


I 


TCS 


NLIE 


TCS 


NLIE 


0.1 


32.21645 


n/a 


18.76030 


n/a 


0.5 


6.671964 


n/a 


4.037716 


n/a 


1.0 


3.662027 


n/a 


2.434501 


2.434431(#) 


1.5 


2.769403 


n/a 


2.027227 


2.027213 


2.0 


2.385451 


n/a 


1.884404 


1.884388 


2.5 


2.190232 


n/a 


1.829459 


1.829456 


3.0 


2.079959 


n/a 


1.809244 


1.809248 


3.5 


2.012596 


n/a 


1.803873 


1.803886 


4.0 


1.968808 


1.968784(#) 


1.804937 


1.804953 


4.5 


1.938895 


1.938889 


1.808633 


1.808658 


5.0 


1.917648 


1.917635 


1.813191 


1.813224 



Table 4.10: The two one-breather states of the Virasoro minimal model Vir{2t,7) perturbed by $(1,3). 
The energy and the volume are normalized to the mass of the kink, which is the soliton of the unrestricted 
sine-Gordon model. Breather #1 has a = while breather #2 corresponds to a = The TCS data 
shown have the predicted bulk energy term subtracted. 



4.10 Conclusions 

In this thesis it is studied how the nonlinear integral equation deduced from the light cone lattice 
model of Q describes the excited states of the sine-Gordon/massive Thirring theory. The most 
important results are summarized as follows: 

1. A derivation of the fundamental NLIE is presented from the light cone lattice which correctly 
takes into account the behaviour of the multivalued complex logarithm function. 



4.10. 



CONCLUSIONS 



87 



2. By examining the infrared limits of the equation it has been shown that (1) it leads to the 
correct two-particle S-matrices for both scattering states and bounded states; (2) it is in 
agreement with the predictions of the TCS method if one chooses the correct quantization 
conditions for the source terms. 

3. By computing the UV conformal weights from the NLIE we have shown that it is consistent 
with the UV spectrum of sG/mTh theory only if we choose the parameter 5 (i.e. the 



quantization rule) as indicated in (4.31 



4. The predictions of the NLIE have been verified by comparing them to results coming from 
the TCS approach (for sG/mTh). 

5. The framework required to deal with minimal models perturbed by 3) is built up. Many 



examples and numerical/analytical checks are given for the ground state in the unitary |32| 
and non unitary cases. 

6. All the conformal dimensions can be reproduced (Kac table) with the convenient choice of 
the twist ( 4.59| ). The particular relation suggested in is not enough to describe the 



whole spectrum. 

7. The IR computations given in section |4.6.4| reproduce correctly the S-matrix of perturbed 
minimal models in the attractive regime (at list in the cases involving only solitons), because 
the attractive second determination drops the twist. 

8. Numerical calculations of concrete examples give a strong evidence for the correctness of the 
energy levels derived from the twisted NLIE for excited states. 

The understanding of "(twisted) sine-Gordon NLIE" and of the finite size behaviour of the con- 
tinuum theory defined from the NLIE ( p. 54 ) is not quite complete. Indeed some open questions 



have not an answer, until now and also further interesting developments can take place: 

1. Is the set of scaling functions provided by the NLIE complete i.e. can we find to every 
sG/mTh or minimal model state a solution of the NLIE describing its finite volume be- 
haviour? This can be called a "counting problem". The main difficulty is that the structure 
of the solutions is highly dependent on the value of the coupling constant - to see that it is 
enough to consider e.g. the appearance of special sources. 

2. The multi-kink states characteristic of the perturbed minimal models (except for the series 
Vir{2,2n + 1)) have been omitted in the analysis so far performed. Although the general 



treatment of the IR spectrum in section 4.6.4 is valid for those states too, a detailed descrip- 



tion is far more complicated than for states which contain only breathers and is left open to 
further studies. 

3. There is also an unresolved technical difhculty, namely that the source configuration of the 
NLIE may change as we vary the volume parameter I. Typically what happens is that while 
the counting function Z is monotonic on the real axis for large volume, this may change as 
we lower the value of I and so-called special sources (and accompanying holes) may appear. 
We do not as yet have any consistent and tractable numerical iteration scheme to handle 
this situation, although the analytic UV calculations and intuitive arguments show that the 
appearance of these terms in the NLIE is consistent with all expectations coming from the 
known properties of perturbed GET. In addition, in the range of I where we can iterate the 
NLIE without difficulty, our numerical results show perfect agreement with TGS. We want 
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to emphasize that these transitions are not physicah the counting function Z and the energy 
of the state are expected to vary analytically with the volume. It is just their description by 
the NLIE with requires a modification of source terms. As it was pointed out also in fs], |9|, 
the whole issue is related to the choice of the branch of the logarithmic term in the NLIE 
(1331). 



The problem pointed out at the previous point is very similar to the behaviour of singularities 
encountered in the study of the analytic continuation of the TBA equation ||l^ and we can 
hope that establishing a closer link between the two approaches can help to clarify the 
situation. From the form of the source terms in the NLIE it seems likely that the excited 
state equations can be obtained by an analytic continuation procedure analogous to the 



one used in TBA [18| to obtain the excited state TBA equations. Certain features of the 
arrangement of the complex roots in the attractive regime and their behaviour at breather 
thresholds also point into this direction. This an interesting question to investigate because 
it can shed light on the organization of the space of states and can lead closer to solving the 
counting problem described above. 



Even if on the lattice the Bethe vectors given in ( 2.37| ) completely describe the Hilbert 



space of the XXZ chain, the form of eigenvectors for continuum energy and momentum is 
completely unknown. This fact reflects the unresolved question of the form of energy and mo- 
mentum eigenvectors for the minkowskian sine-Gordon theory (the Faddeev-Zamolodchikov 
algebra is only a phenomenological picture) . Also the determination of correlation functions 



is in an early stage, but in recent publications |46| explicit expressions for the field ^e*"*^) 
and its descendent are given. 

6. The extension of light-cone approach and NLIE to other QFT models is in progress. An 
interesting extension is the description of finite volume spectrum of Ziber-Mikhailov-Shabat 
model (imaginary Bullogh-Dodd). A first suggestion in this direction, even if is obtained in 



a completely different approach, recently appeared in the |^, for the ground state. 



The second and fourth points are important because their investigation may lead closer to un- 
derstanding the relation between the TBA and the NLIE approaches. It is quite likely that 
establishing a connection between the two methods would facilitate the development of both and 
may point to some common underlying structure. 

As a final remark, I want to speak about an application of the NLIE to a chemical compound. 



i.e. "copper benzoate" (Cu(C6D5COO)2 • 3D2O). Its specific heat has been computed in [|45| 
using the vacuum untwisted (a = 0) NLIE, where the finite size L is the inverse temperature 
L = {kBT)~^. This corresponds to do equilibrium thermodynamics of sine-Gordon model. 

A curiosity: in that contest the NLIE is called Thermal Bethe Ansatz, emphasizing the well 
known relation between statistical mechanics and QFT. Exactly in a similar contest (Heisenberg 
ferromagnet model) was born the Bethe Ansatz 



Appendix A 

Fourier transformation: some 
conventions 



Given a bounded continuous function f{x) defined on the whole real axis, its Fourier transform 
is defined as (all the integrals must be taken on the whole real axis): 



1 

-ikx 



dxe P+'^f{x) (A.l) 
and the inverse Fourier transform is 

1 

The convolution of two functions, defined as follows, can be expressed in terms of the Fourier 
transforms of the two functions: 

1 

/I r .11, ike 
dx f{e - x) g{x) = —^J-e P + 1 fikygik) (A.3) 

The following expression holds between the Fourier transforms of a function and its derivative: 
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Appendix B 

The function (/)(A, iy) 



In this appendix all the important properties of the function (f) will be clarified. Let 

sinh (iiTU + t?) 

=i\og , = V i9eC (B.l) 

sinh (ivrzv — {)) 

and continuity is required around the real axis. The interest is for = 1/2, 1 and p > 0. The 
requirement of oddity implies, for real values of that (p{0, u) = 0. Oddity and continuity can be 
implemented only if the fundamental determination (FD) of the logarithm is assumed, in a strip 
containing the real axis: 

— TT < 9m(logp£) w) < TT. 

This choice implies that (p is real on the real axis. The argument of the log has poles and zeros. 
They are essential singularities for the cp function. Their position is given by: 

9mi? = ibvr (A;(p+ 1) - I/) , k e Z. ^ 
The fundamental strip around the real axis (FD) is bounded by the first encountered singularity: 
■i? G M X ]-min{z^7r, it (p + 1 - i^)} ,min{z^7r, n (p + 1 - (B.3) 



To extend the definition of ([B.lJ ) to the whole plane it is necessary to give a prescription for the 
position of the cuts (because of the logarithm, in a small closed trip around one of the singularities 



, the value of (p changes by 2tt). This is done as in figure [B.l 



The derivative of (p is the single value meromorphic function: 

2 sin ■ 



cosh cos 

p+l p+l 

This shows that, on the real axis, 4> is monotonically increasing if < — -j-y < tt/2 and decreasing 

if 7r/2 < < TT. The asymptotic values, for in the fundamental strip, are: 

p+l 

lim </,(^,z.)=±vr(l-^). (B.5) 

'0^±OO P+l 
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1 


/ ttT n+1 - 

-\--- 


-2n 


-2k 


--4 


- iKV 

4-- 


-2k 

-/7r(p+l-v) 


-2k 



1 




+271 






/7r(p+l-v) 


\-- 

+2k 


t 


-/7r(p+l-v) 


+2k 






- iKV 


1--- 

+271 



(a) 



(b) 



Figure B.l: Positions of singularities and cuts for the function 4>. (a) in the case — -—^ < 1/2; (b) in the 

case > 1/2. The little arrows indicate the jump in the values of the function. 

p + 1 



Out of the fundamental strip, the prescription indicated in figure |B.lJ must be used. 

The Fourier transformation of (j)' can be obtained from the definition (A.l), using the theorem 
of residues: 



sinh ■ 



(l)'{k,u) = 27r- 



TT f p + 1 — 21^ 
P+l 



k 



sinh —k 



(B.6) 



This formula is the Fourier transform of 4>'{'d, v) only for d in the fundamental strip ([B.3|) , because 



out of this strip new singularities of (B.4) appear in the computation of Fourier transform. 



Appendix C 

A lemma for UV computations 



In the following lemma has been proved. Assume that f{x) satisfies the non-linear integral 
equation 

/dy 
— G{x — y — ie) logfl + fix + ie)) 

where ^{x) is real on the real axis and G{x) = G{—x) is real too, with bounded integral and 
peaked around the origin. From this equation follows that f{x) has unit modulus on the real axis. 
To avoid crossing of the branch cut of the logarithm, assume that 

if /(x + ie)GM then f{x + ie)>-l. (C.l) 

Then the following expression holds: 

/dx 
— if'ix + ie) logfl + fix + ie)) = 

= -2 3f?e j ^log(l + u)-^ j dxG{x) {F\+oo) - F\-oo)) 



(C.2) 



where the curve T is any path in the complex plane that goes from /_ = /(— oo + ie) to /+ = 
/(-|-oo + ie) avoiding the branch cut, i.e. respecting the condition (CJ.) and 

F{x) = 2 9m log(l + f{x + ie)). 

In the 6- vertex case all the terms can be made manifest by the identification of G with that defined 
in ( ^^211) , f{x) = (-l)'5e*^±W and ip{x) with (14251) . 

Consider the kink "+". The integral on the path T goes from the points 

/_ = f{-oo + ie) = (-l)'5e*^+(-°°) = e*2+(-°°) 
= /(+CX) + ie) = (_i)'5e*^+(+~+»^) = 

(the first one can be computed with (|4.17D ). Observe that |/_| = 1 then the path F can be 
composed by an unit radius arc from /_ to the point 1 and a segment from 1 to 0. On the arc the 
integration variable has unit modulus then it is convenient to do the change of variables u = e*". 
This simple computation gives 



2SRe f —log{l + u) = -Qli-oo). 

J arc ^ 
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The integration on the segment gives a well known dilogarithmic expression: 

^2 



_25?e r^log(l + n)=2 ^ log(l + u) = ^ 
Ji u Jq u 6 

The quantities F(iboo), from their definition, give: 

F(+oo) = 

F(-oo) = 2 9m log(l + (_i)5ei^+(-oo)) ^ g^(_oo). 

Remembering now that the integral of G appearing in ( p. 21 ) is expressed in ( p. 33 ) yields (the 
computation for the "-" kink is completely analogous): 



2Qm J y(^±(x + ir,) log (l + (_i)'5ei^±(x'+i'7) 



(C.3) 

6^4 p 
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